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Abstract.  Plants  and  animals  that  inhabit  the  intertidal  zone  of  wave-swept 
shores  are  generally  small  relative  to  terrestrial  or  subtidal  organisms. 

Various  biological  mechanisms  have  been  proposed  to  account  for  this  observation 
(competition,  size-specific  predation,  food-limitation,  etc.).  However,  these 
biological  mechanisms  are  constrained  to  operate  within  the  mechanical 
limitations  imposed  by  the  physical  environment,  and  these  mechanical 
limitations  have  never  been  thoroughly  explored.  It  is  thus  possible  that  the 
observed  limits  to  size  in  wave-swept  organisms  are  due  solely  or  in  part  to 
mechanical,  rather  than  biological,  factors.  This  possibility  is  examined. 

The  total  force  imposed  on  an  organism  in  flowing  water  is  due  both  to  the 
water's  velocity  and  to  its  acceleration.  The  force  due  to  velocity  (a  combined 
effect  of  drag  and  lift)  increases  in  strict  proportion  to  the  organism's 
structural  strength  as  the  organism  increases  in  size,  and  therefore  cannot  act 
as  a  mechanical  limit  to  size.  In  contrast,  the  force  due  to  the  water’s 
acceleration  (the  acceleration  reaction)  increases  faster  than  the  organism's 
structural  strength  as  the  organism  grows,  and  thus  forms  a  potential  mechanical 
limit  to  its  size.  This  fact  is  incorporated  into  a  model  that  predicts  the 
probability  that  an  organism  will  be  destroyed  (by  breakage  or  dislodgement)  as 
a  function  of  five  parameters  that  may  be  empirically  measured: 

1.  The  organism's  size, 

2.  The  organism's  structural  strength, 

3.  The  maximum  water  acceleration  in  each  wave, 

4.  The  maximum  water  velocity  at  the  time  of  maximum  acceleration  in  each 
wave,  and 

5*  The  probability  of  encountering  waves  with  given  flow  parameters. 

The  model  is  tested  using  a  variety  of  organisms.  For  each,  parameters  1-4 
are  measured  or  calculated; -the  probability  of  destruction,  and  the  size- 
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specific  increment  in  this  probability,  are  then  predicted.  In  the  cases  of  the 
limpets,  Collisella  pelta  and  Notoacmaea  scutum,  the  urchin,  Stronisylocentrotus 
purpuratus,  the  mussel,  Mytilus  calif orni anus  (when  solitary),  and  the 
hydrocoral,  Millepora  complanata,  the  probability  of  destruction  and  the  size- 
specific  increase  in  the  risk  of  destruction  are  both  substantial.  It  is 
conjectured  that  the  size  of  individuals  of  these  species  may  be  limited  as  a 
result  of  mechanical  factors,  thou^i  the  case  of  complanata  is  complicated  by 
the  possibility  that  breakage  may  act  as  a  dispersal  mechanism.  In  other  cases 
(the  snails,  Thais  canaliculata,  Thais  emarginata,  and  Littorina  scutulata;  the 
barnacle,  Semibalanus  cariosus)  the  size-specific  increment  in  the  risk  of 
destruction  is  small  and  the  size  limits  imposed  on  these  organisms  are 
conjectured  to  be  due  to  biological  factors. 

Aside  from  providing  evidence  regarding  the  mechanical  limits  to  size  of 
wave-swept  organisms,  the  model  constructed  here  provides  an  experimental 
approach  to  examining  many  potential  effects  of  environmental  stress  caused  try 
flowing  water.  For  example,  these  methods  may  be  applied  to  studies  of: 

1.  Life  history  parameters  (e.g.  size  at  first  reproduction,  age  at  first 
reproduction,  timing  of  reproductive  cycles,  the  length  of  the  possible 
reproductive  lifetime), 

2.  The  effects  of  gregarious  settlement  on  the  flow  encountered, 

3.  The  physical  basis  for  patterns  of  disturbance, 

4.  The  optimum  (as  opposed  to  the  maximum)  size  of  or^inisms,  and 

5.  The  energetic  cost  of  maintaining  a  skeleton  with  an  appropriate  safety 
factor. 

A  definitive  answer  regarding  the  possibility  of  mechanical  limits  to  size 
depends  both  upon  an  accurate  measurement  of  encountering  a  wave  of  specific 
flow  parameters  and  upon  factors  thau  are  external  to  the  model  considered  (e.g. 
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life  history  parameters).  Further,  due  to  their  ability  to  move  with  the  flow, 
organisms  that  are  sufficiently  flexible  can  escape  the  size  limits  imposed  on 
more  rigid  organisms.  Thus,  some  macroalgae  can  attain  large  sizes  (2-3  m  in 
maximum  dimension).  The  precise  role  of  these  factors  awaits  further  research. 

Key  Words,  mechanical  limits,  size  limits,  wave  forces,  intertidal  organisms, 
corals,  acceleration  reaction,  added  mass,  size-specific  mortality,  exposure, 
disturbance. 

10  Year  Index  Bitries.  Mechanical  limits  to  size  of  wave-swept  organisms, 
quantitative  measurement  of  exposure,  size-specific  incremental  risk  of 
dislodgament,  size-specific  mortality,  size-specific  force  due  to  water 
acceleration,  probability  of  dislodgment,  contributions  of  biomechanics  to 
ecology. 
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INTRODUCTION 


Consequences  of,  and  constraints  on,  the  size  of  living  things  have  long 
intrigued  observers  of  nature  (g.g.,  Galilei  1638,  Thompson  1917,  Haldane  1928). 
Various  patterns  have  been  noted  in  the  ecological  roles  of  organisms  of 
different  sizes  (Hutchinson  and  MacArthur  I960,  VanValen  1975,  Kerr  1974,  Banse 
and  Mosher  1980,  Silvert  and  Platt  1981,  but  see  Frost  1980)  and  in  the  size 
distribution  of  organisms  through  evolutionary  time  (g.g.,  Bonner  1968,  Gould 
1966,  Stanley  1973).  Basic  "rules"-  for  the  dependence  on  body  size  of  the  rates 
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at  which  organisms  function  (e.g.  respire,  capture  food,  grow,  or  reproduce) 
have  been  noted,  and  various  explanations  have  been  suggested  (see  Schmidt- 
Nielsen  1 974,  McMahon  1 975a,  1 980,  Banse  and  Mosher  1 980,  Gray  1 981 ,  Platt  and 
Silvert  1981,  and  references  cited  therein).  Organisms  of  grossly  different 
size  often  live  in  different  physical  realms  (g.g.,  Haldane  1928,  Went  1968, 

Horn  et_  ah_  1982).  From  the  many  important  questions  related  to  the  size  of 
organisms  we  have  chosen  to  examine  the  mechanical  factors  which  may  define  the 
upper  boundary  to  size. 

Upper  Limits  to  Body  Size 

Factors  imposing  physical  upper  limits  to  size  have  received  considerable 
attention.  Since  the  volume  of  an  organism  is  proportional  to  [length]^,  whereas 
the  surface  area  for  exchange  of  materials  as  well  as  the  cross-sectional  area 
of  skeletal  elements  (and  hence  their  strength)  are  both  proportional  to 
[length]^,  the  shapes  of  many  large  organisms  are  different  from  small  ones  in 


ways  that  permit  area  to  keep  pace  with  volume  (e.g.,  Haldane  1928,  Thompson 
1917,  Gould  1966,  Stanley  1973,  Schmidt-Nielsen  1974).  Such  allometric  changes 
in  morphology  with  increase  in  body  size  eventually  lead  to  structural 
absurdities  (Gould  1966):  at  large  body  sizes  skeletons  would  deflect  too 
much,  break,  buckle,  or  be  too  bulky  to  be  moved  (e.g.,  Galilei  1638,  Haldane 
1928,  Gould  1966,  Currey  1970,  MacMahon  1973).  The  mechanical  practicality  or 
Ibe  metabolic  cost  of  various  types  of  locomotion  (running,  swimming,  flying, 
burrowing)  may  set  physical  limits  to  the  sizes  of  various  organisms  (e.g., 

Gould  1966,  Bonner  1968,  Maynard-Smith  1968,  Alexander  1971,  Stanley  1973, 
Schmidt-Nielsen  1974,  for  examples  see  Pedley  1977).  And,  as  we  discuss  later, 
many  other  factors  have  also  been  proposed  to  limit  the  size  of  organisms. 

Here  we  examine  size  limits  of  organisms  that  live  in  wave-swept 
environments,  in  particular  those  organisms  that  inhabit  the  intertidal  zone  of 
rocky  shores  and  coral  reefs.  Wave-swept  shores  are  undoubtedly  a  physically 
harsh  environment,  and,  as  such,  should  form  an  ideal  testing  ground  for  an 
examination  of  the  role  played  by  the  physical  environment  in  limiting  size,  and 
thereby,  for  the  role  played  by  this  one  aspect  of  the  environment  in  the 
ecology  of  these  organisms. 


Wave-Swept  Organisms  Are  Not  large 

One  simple  observation  forms  the  basis  for  this  examination:  Benthic 
organisms  that  live  near  the  water's  surface  on  wave-swept  shores  do  not  attain 
large  size.  Not  only  are  the  largest  of  wave-swept  organisms  (seaweeds  and 
corals,  which  rarely  exceed  2  to  3  m  in  their  greatest  dimension)  much  smaller 
than  the  largest  organisms  on  the  face  of  the  earth  (whales;  Sequoia  trees 
(Alexander  1971)),  but  they  are  also  smaller  than  the  largest  benthic  marine 


organisms  living  in  more  protected  habitats.  This  trend  has  been  observed,  for 
example,  on  rocky  temperate  shores  (lewis  1968).  Mussels  are  larger  in  those 
intertidal  areas  where  wave  action  is  lower  (Harger  1970,1972),  or  subtidally 
(Baine  1976a, 1976b)  where  the  intensity  of  water  motion  is  less.  The  same  trend 
has  been  noted  for  starfish  (Paine  1976a).  Many  algae  are  stunted  in  hi^i  wave 
exposure  (Schwenke  1971,  Connell  1972),  although  there  are  exceptions  to  this 
trend  (Schwenke  1971).  Colonial  animals  such  as  gorgonians,  hydroids,  corals  and 
sponges  produce  taller  colonies  in  deeper  water  where  wave  action  is  attenuated 
(Reidl  1971).  Similar  trends  in  maximum  body  size  have  been  noted  in  tropical 
coral  reef  systems.  On  those  reefs  or  portions  of  reefs  most  battered  by  waves, 
corals  are  smaller  or  have  lower  profiles  than  on  more  protected  sections  of 
reefs  (Stoddart  1969,  Glynn  1973,  Endean  1976,  Randall  and  Eldredge  1977,  Smith 
and  Harrison  1977,  Adey  1978,  Highsmith  1981).  Vithin-species  comparisons  show 
the  same  trend  for  stony  corals,  gorgonians,  and  algae  (Grassle  1973,  Vosburgh 
1977,1982). 

Why  should  wave-swept  organisms  be  limited  to  these  small  sizes?  There  are 
many  biological  factors  which  may  be  responsible  for  the  small  sizes  of  these 
organisms  (food  limitation,  size-specific  predation,  etc.),  and  these  factors 
have  been  addressed  by  many  of  the  studies  of  intertidal  ecology  cited  above. 
However,  biological  interactions  must  operate  within  mechanical  limits  set  by 
the  physical  environment.  These  limits  have  never  been  defined  for  wave-swept 
organisms,  and  it  is  thus  quite  possible  that  apparent  biological  limits  are 
artefacts.  Consequently,  in  this  study  we  postpone  until  the  Discussion  the 
consideration  of  biological  interactions,  and  examine  instead  the  hypothesis 
that  there  are  purely  mechanical  limits  to  the  size  to  which  organisms  may  grow 
in  wave-swept  environments.  There  is  considerable  evidence  which  supports  this 


hypothesis,  showing  that  large  organisms  are  indeed  more  likely  to  be  ripped  off 
hard  substrata  by  waves  than  are  small  ones  (e.g.  corals:  Graus  et_  aL_  1 977, 
Randall  and  Eldredge  1977,  Adey  1978,  Chamberlain  1978,  Higjismith  1980,  1981, 
1982,  mussels:  Harger  1970,  Paine  and  Levin  1981,  algae:  Black  1976,  Santelices 
et  al.  1 980,  and  other  rocky  shore  organisms,  see  Connell  1 972).  Once  this 
hypothesis  has  been  evaluated,  the  role  played  by  biological  interactions  can  be 
examined  more  objectively. 

The  Role  of  Acceleration 

The  basis  of  our  argument  is  a  consequence  of  the  unsteady,  accelerational 
nature  of  the  flow  in  wave-swept  habitats.  Most  past  studies  on  the  forces 
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exerted  on  wave-swept  organisms  have  considered  only  the  hydrodynamic  forces  due 
to  steady,  nonaccelerating  water  flow  (e.g.  Jones  and  Demetropoulos  1968,  Branch 
and  Marsh  1978,  for  exceptions  see  Carstens  1968,  Koehl  1977a).  This  force  is 
the  drag  force, 

ZdocA^IJ2  0) 

where  A^  is  the  area  of  the  organism  exposed  to  flow,  projected  in  the  direction 
of  flow,  and  U  is  water  velocity*  This  force  is  applied  to,  and  resisted  by, 
some  critical  structure  of  the  organism  having  an  area  AqJ  for  example  the  basal 
attachment  area  of  an  acorn  barnacle.  Thus  the  stress  (force/area,  cr)  placed  on 
the  organism's  structural  material  is: 

&  ~  gd/A^tt  (Ap/ArJ  U2 


(2) 


?or  an  organism  that  maintains  a  set  shape  as  it  grows  (i-e.  grows 
isometrically)  A^/Aq  is  constant  and 

(5  =  B  U2  (3) 

where  B  is  the  appropriate  constant  of  proportionality.  It  is  evident  from  eq.  3 
that  for  steady  flow  the  stress  placed  on  an  organism  is  independent  of  its 
size.  Certainly,  the  larger  the  organism,  the  larger  the  force  to  which  it  is 
exposed,  hut  the  area  over  which  this  force  is  applied  is  proportionally  larger 
as  well,  so  the  stress  is  constant.  Thus  the  hydrodynamics  of  steady  flow  give 
no  clue  to  size  limitation  in  isometrically  growing  organisms  as  long  as  there 
is  no  significant  spatial  variation  in  U  such  that  bigger  organisms  encounter 

faster  flow. 

\ 

However,  water  flow  in  waves  is  not  steady.  Water  alternately  moves  towards 
and  away  from  shore  and  is  thus  periodically  accelerated.  When  waves  break, 
large  accelerations  accompany  the  turbulent  eddies  that  are  formed.  Objects  in 
accelerating  flow  experience  a  force,  £3,  the  acceleration  reaction,  in  addition 
to  the  forces  caused  by  the  instantaneous  water  velocity.  A  more  complete 
(though  less  than  exhaustive)  expression  for  the  total  force  imposed  on  an 
organism  is 

*df  =  (?d  +  4)  cc  (A^  U2)  +  (V  dU/dt)  (4) 

where  is  the  total  force  in  the  direction  of  flow  and  V  is  the  volume  of 
water  displaced  by  the  organism  (Batchelor  1967,  Daniel  1982).  If  we  divide  this 
equation  by  the  critical  area  over  which  the  force  is  applied  and  assume  that 
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the  organism's  volume  =  L  where  L  is  thus  defined  as  a  characteristic  length, 
the  stress  on  the  organism  is: 

cr=  ?df/Ac  =  B  U2  +  B*  L  dU/dt  (5) 

where  B  and  B'  are  the  appropriate  constants  of  proportionality.  Thus,  in  an 
accelerating  fluid  the  stress  experienced  by  an  organism  is  indeed  a  function  of 
the  organism's  characteristic  length:  the  larger  the  organism,  the  greater  the 
stress.  For  any  given  accleraticnal  flow  if  L  exceeds  a  certain  value,  exceeds 

the  stress  that  the  organism  can  sustain  (the  breaking  stress)  and  the  ganism 

dislodges  or  breaks.  In  this  manner  accelerating  flow  acts  to  set  a  pc  tial 
mechanical  limit  to  size. 

In  order  to  examine  this  possibility  we: 

1.  Create  a  model  for  the  stresses  placed  on  wave-swept  organisms, 

2.  Estimate  the  flow  regime  present  on  wave-swept  shores, 

3-  Measure  the  necessary  hydrodynamic  coefficients  and  constants  for  actual 
wave-swept  organisms, 

4.  Measure  the  distribution  of  adhesive  tenacity  for  various  wave-swept 
species,  and 

5.  Apply  the  principles  of  the  model  to  these  species. 


MATERIAIS,  METHODS,  AND  INTERMEDIATE  RESULTS 


The  exposed  rocky  shoreline  of  Washington  State  and  a  coral  reef  in  the 
western  Caribbean  were  chosen  as  model  systems  in  which  to  test  this  theory  of 
size  limitation.  The  majority  of  the  empirical  measurements  made  were  carried 
out  at  various  sites  on  Tatoosh  Island,  Washington  (124°44'W,  4802311).  The 
principal  site  is  a  rock  face  inclined  approximately  35°  from  horizontal  and 
fully  exposed  to  to  the  prevailing  southwesterly  swells.  The  bottom  slopes 
steeply  down  from  the  site  reaching  a  depth  of  10m  (below  MLLW)  approximately 
30m  offshore.  This  depth  is  maintained  for  another  70m  out  from  the  site  before 
the  bottom  again  slopes  steeply  downward.  Depths  of  80-1 00m  are  reached  within 
500m  of  the  site.  Other  measurements  were  conducted  in  various  surge  channels 
on  Tatoosh.  Further  intertidal  measurements  were  made  at  Mukkaw  Bay 
(approximately  three  kilometers  south  of  Tatoosh  Island),  and  at  Shi-Shi  Beach 
(approximately  8  kilometers  south  of  Tatoosh  Island),  in  both  cases  on  rocks 
protected  from  the  full  force  of  the  prevailing  swells  by  the  gently  sloping 
bottom  offshore.  Measurements  of  wave  height,  coral  morphology  and  coral 
distribution  were  made  at  the  Galeta  Point  Marine  laboratory  of  the  Smithsonian 
Tropical  Research  Institute  near  Colon,  Republic  of  Panama  (78°43,W,9°34'N).  As 
with  any  intertidal  site,  the  fine-scale  topography  of  the  sites  used  in  this 
study  is  so  complex  that  a  detailed  description  is  of  little  use  for  fluid 
mechanical  purposes. 

Estimation  of  Water  Velocity  and  Acceleration 


December  1979  and  November  1980-  A  wave-force  telemetry  system  was  placed  in  the 
rock  surface  as  described  by  Denny  (1982)  and  the  magnitude  and  direction  of 
tydrodynamic  forces  exerted  on  cast  plastic  models  of  a  limpet  (Collisella 
pelta)  and  on  spheres  (1  cm  diameter)  were  measured  for  each  wave  during  the 
course  of  several  tides.  Approximately  3  x  104  force  events  were  recorded.  The 
data  for  this  study  were  taken  from  the  portion  of  the  record  where  the  most 
violent  waves  occurred  (November  19-20,  I960).  During  this  period  swell  heights 
at  breaking  were  visually  estimated  to  vary  from  2  to  4  m.  ?or  each  of  30 
sequential  large  waves  the  time  vs.  force  record  was  used  to  calculate  time  vs. 
velocity  and  time  vs.  acceleration  records  by  numerically  solving  a 
rearrangement  of  eq  4: 


an/dt  =  (2af-[p/2]  ^  u2)/!^  v] 


(6) 


using  a  fourth  order  Runge-Kutta  method  with  the  boundary  condition  that  all 
force  at  zero  time  is  due  to  the  acceleration  reaction.  ?or  the  derivation  of 
this  equation  see  eqs.  1 7-1 9  below.  The  volume,  V,  and  projected  area,  A^  of  the 
limpet  replica  used  in  these  measurements  were  determined;  and  the  drag 
coefficient,  £d»  and  inertia  coefficient,  C^,  were  empirically  determined  as 
discussed  below.  The  density  of  seawater,  p,  is  assumed  to  be  constant  at  1025 
kg/rn^. 


Estimation  of  Extreme  Wave  Conditions 


accurately  "the”,  extreme  value  for  a  randomly  occurring  phenomenon  such  as  a 
breaking  wave.  ?or  example,  during  the  course  of  this  study  low  tides  suitable 
for  installing  the  telemetry  system  did  not  coincide  with  any  severe  storms,  and 
the  records  obtained  are  thus  of  waves  of  less  than  maximum  size.  On  the  other 
hand,  it  would  have  been  hazardous  and  probably  technically  impossible  to 
install  the  device  during  more  extreme  environmental  conditions.  The  maximum 
values  of  velocity  and  acceleration,  calculated  as  described  above,  are 
approximately  14  m/s  and  400  m/s^.  This  maximum  calculated  acceleration  occurred 
at  a  calculated  velocity  of  1.7  m/s,  and  thus,  assuming  that  the  water  started 
with  zero  velocity,  the  force  accelerating  the  water  was  applied  for  at  most 
about  4  ms.  The  response  time  of  the  transducer  used  is  8  ms  (Denny  1982),  so  it 
is  very  likely  that  the  transducer's  deflection  lagged  behind  the  force  applied 
and  that  the  forces  measured  here  are  underestimates.  Consequently  the 
velocities  and  accelerations  calculated  will  be  underestimates  of  the  actual 
values  for  the  relatively  mild  conditions  encountered  and  thus  may  considerably 
underestimate  the  conditions  accompanying  storm  waves-  Rather  than  guessing  at 
one  specific  maximum  velocity  and  acceleration  we  use  the  values  calculated  here 
as  a  rough  guide  to  the  lower  limit  of  what  the  actual  maxima  are,  and  present 
data  for  velocities  from  0  to  20  m/s  and  accelerations  from  0  to  2000  m/s^.  At 
present  we  have  no  independent  evidence  to  confirm  this  range  of  accelerations, 
however  there  is  confirming  evidence  for  the  range  of  velocities.  Naval  records 
of  wave  heists  (U.S.  Navy  1973)  show  that  for  deepwater  waves  recorded  off 
southwestern  Vancouver  Island  (near  the  Tatoosh  Island  site)  between  1949  and 
1973  the  maximum  wave  height  was  7.75  -  9-75  m  on  0.1$  of  the  days  when  heights 
were  recorded,  and  no  waves  higher  than  this  were  recorded.  As  will  be  shown 

below  (eq.  14),  the  estimated  maximum  water  velocity  at  the  moment  the  wave 
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breaks  for  a  9.75  m  wave  is  approximately  14  m/s,  although  this  may  increase  as 
the  wave  crest  plunges  downwards  (Carstens  1968).  Further,  the  shoaling  of  waves 
as  they  near  the  shore  may  increase  their  height  above  that  observed  in  deep 
water.  In  ligit  of  these  factors  the  velocity  range  used  here  seems  justified. 

Measurement  of  Drag  and  Lift  Coefficients 

Drag  coefficients  were  measured  for  two  standard  shapes  (sphere, 
hemisphere)  and  three  species  of  locally  abundant  and  characteristic  shelled 
animals  (an  acorn  barnacle,  Semibalanus  cariosus,  a  limpet,  Collisella  pelta, 
and  a  snail,  Thais  (Nucella)  canaliculata).  Lift  coefficients  were  measured  for 
S.  cariosus  and  pelta.  The  object  to  be  tested  was  attached  to  the  mounting 
bolt  of  a  force  transducer  sensitive  either  to  the  force  in  the  direction  of 
water  movement,  F^f,  or  the  force  perpendicular  to  water  movement,  i.e.  lift, 
jjj..  The  force  transducer  was  mounted  on  the  top  wall  of  a  recirculating  water 
tank  (similar  to  those  described  by  Vogel  and  LaBarbera  1978),  the  mounting 
bolt  extending  througi  a  hole  in  the  wall  such  that  the  object  was  held  within  1 
-  2  mm  of  the  wall.  Each  object  was  then  subjected  to  steady  water  velocities  up 
to  approximately  4  m/s.  Assuming  the  flow  to  be  fully  turbulent,  the  boundary 
layer  in  the  vicinity  of  the  model  is  calculated  to  be  <  5  mm  (see  equation  15). 
Mainstream  velocities  in  the  test  section  of  the  flow  tank  were  measured  with  an 
electromagnetic  flowmeter  (Cushing  Velmeter  600P).  The  projected  area  of  each 
object  (in  the  direction  of  the  force  being  measured)  was  determined  by  tracing 
an  appropriately  projected  picture  of  the  object  onto  a  piece  of  paper.  The 
resulting  area  was  cut  out,  weighed,  and  compared  to  the  we i git  of  standard 
areas  to  arrive  at  the  final  estimate  of  projected  area  For  each  object  the 


force  was  measured  at  several  water  velocities  and  the  lift  and  drag 
coefficients  were  calculated  using  the  following  equations  (Vogel  1 981 ): 

£d  =  2?d/(f^!I2)  (V) 

=  2FjJ(ekR  U2)  (8) 

In  general  both  and  Cj_  vaiT  with  object  size,  the  physical 
characteristics  of  the  fluid,  and  the  fluid  velocity.  In  order  to  compare 
coefficients  between  objects  of  different  sizes  and  Cj  were  plotted  as  a 
function  of  Reynolds  number,  Re,  a  dimensionless  value  which  is  a  measure  of  the 

i 

relative  contribution  of  inertial  and  viscous  processes  in  determining  flow 
patterns  around  objects: 

Re  =  U  L/  v  (9) 

where  U  is  the  velocity  of  the  fluid  relative  to  the  organism,  L  is  a 
characteristic  length  which  we  take  to  be  the  diameter  of  the  organism,  and 
is  the  fluid's  kinematic  viscosity  (1  x  .10“^  m  ^/s)  (?igure  1).  The  values  used 
here  are  those  obtained  for  Reynolds  numbers  around  10^  (equivalent  to  a 
velocity  of  10  m/s  for  an  object  1  cm  long),  a  value  generally  approriate  for 
the  conditions  encountered  during  breaking  waves.  A  variety  of  water  velocities 
and  sizes  of  objects  are  used  in  this  study,  corresponding  to  Re's  both  hinder 
and  lower  than  the  standard  Re  used  to  estimate  and  Cj_.  Thus  the  values  may 
slightly  over-  or  underestimate  the  actual  and  va^-ues  i°  given 
situation.  This  dependence  of  and  C-^  on  Re  (and  thereby  on  size)  is  likely  to 


be  most  evident  for  simply  shaped,  smooth  organisms  such  as  the  limpets  measured 
here.  As  Vogel  (1 981 )  observes,  the  variation  of  with  Fte  decreases  as  the 
irregularity  of  the  object  increases;  thus  for  an  organism  such  as  a  sea  urchin, 
Cjj  may  vary  little  with  Re.  The  C<i  of  plate-like  organisms  (such  as  the  coral 
species  discussed  later)  varies  only  insignificantly  (<5$)  over  the  Re  range  of 

10^-10^  (Hoerner  1965).  The  effect  of  variations  of  and  Cq  is  examined 
more  closely  in  the  Discussion  and  in  Appendix  2.  The  and  values  used  here 
were  measured  in  steady  flow.  If,  instead,  the  flow  oscillates,  and  may 
vary  as  a  function  of  another  dimensionless  parameter,  the  period  parameter,  K 
(Sarpkaya  and  Isaacson  1 981 ). 

K  =  T^TA  (10) 

where  Um  is  the  maximum  velocity*  T  is  the  period  of  oscillation,  and  L  is  a 
characteristic  length  (usually  the  organism's  diameter).  No  attempt  was  made  to 
measure  changes  in  and  as  a  function  of  K. 

Drag  and  lift  coefficients  for  objects  for  which  empirical  measurements 
were  not  made  have  been  estimated  either  using  the  values  presented  by 
Hoerner(1 965)  or  by  assuming  values  equal  to  those  of  an  object  of  similar 
shape.  ?or  example,  values  for  the  acorn  barnacle  Balanus  glandula  are  assumed 
to  be  equal  to  those  of  cariosus  for  which  measurments  were  obtained  (Table 
1). 

Location  of  Center  of  Pressure 

The  lift  imposed  on  an  organism  is  the  result  of  a  pressure  differential 
across  the  organism.  This  distributed  lift  force  behaves  as  if  the  total  lift 
were  acting  through  a  single  point,  the  center  of  pressure.  The  center  of 


pressure  was  determined  as  follows:  Three  hollow  cones  (basal  diameter  11  cm, 
height/radius  =  0.5,  1.0,  2.0  respectively)  were  turned  from  acrylic  plastic.  In 
each,  five  ports  (approximately  1  mm  diameter)  were  drilled  in  a  line  running 
from  the  basal  edge  up  to  and  including  the  apex  (see  Figure  2).  Each  cone  was 
attached  to  the  top  wall  of  the  flow  tank  and  exposed  to  a  steady  velocity  of 
0.63  m/s  (Re  =  6.9*1 0^).  The  base  of  the  cone  was  not  sealed  to  the  wall, 
allowing  hydrostatic  pressure  to  be  transmitted  to  the  fluid  beneath  the  cone. 
The  hollow  interior  of  the  cone  was  loosely  packed  with  glass  wool  to  inhibit 
gross  flow.  The  pressure  difference  between  each  port  and  a  reference  port 
located  beneath  the  cone  was  measured  manometrically  using  CCI4  =  1594 
kg/m^)  as  the  manometer  fluid.  Readings  were  taken  with  the  row  of  ports  at 
angles  to  the  water  flow  varying  every  3CP  between  0°and  180°  (CP  =  directly 
into  the  flow).  Pressure  readings  were  then  summed  to  estimate  the  location  of 
the  center  of  pressure,  C^: 


wo  c=0  f  ' 
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The  definition  of  terms  is  shown  in  Figure  2. 

The  center  of  pressure  was  found  to  lie  along  the  midline  and  downstream  of 
the  apex  by  distances  varying  between  0.111  Lj.  (for  the  cone  with  height/radius 
=  2)  and  0.29  Lp  (for  the  cone  with  height/radius  =  0.5),  where  Lp  is  the  cone 
radius.  The  average  location  of  the  center  of  pressure  was  0.24  L*.  downstream  of 
the  apex,  and  this  value  is  used  in  the  calculations  here. 
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Measurement  of  Added  Mass  Coefficient 

The  relationship  between  water  acceleration  and  the  force  encountered  by  an 
organism  depends  in  part  on  the  shape  of  the  organism.  This  effect  of  shape  is 
accounted  for  by  use  of  an  added  mass  coeficient, .  This  coefficient  was 
determined  as  follows:  Negative  molds  were  prepared  by  embedding  shells  of  S. 
cariosus,  C.  pelta,  and  T^_  emarginata  in  silicone  sealant.  Following  removal 
of  the  shells,  the  molds  were  injected  with  urethane  foam  (Insta-?oam  Products, 
Inc.)  to  provide  rigid,  light  models.  The  models  were  coated  with  spray  paint  to 
prevent  waterlogging  and  the  mass  of  each  model  was  measured  on  an  electric 
balance  to  +  0.001  g.  The  mass  of  water  each  model  displaced  when  immersed  was 
determined  by  measuring  the  bouyant  force  exerted  by  submerged  models.  The  mass 
of  displaced  fluid  was  measured  to  +  0.05  g. 

Models  were  attached  to  the  tip  of  a  thin  aluminum  rod  that  was  mounted  on 
a  force  platform  similar  to  that  of  Denny  (1982). 

Acceleration  of  the  model  was  accomplished  by  suspending  the  force  platform 
from  the  ceiling  with  wires  approximately  5  m  long  and  allowing  the  platform  to 
swing  like  a  pendulum.  The  model  and  rod  were  at  all  times  immersed  in  a  long 
trough  of  water.  Acceleration  of  the  model  was  measured  using  an  accelerometer 
(Ehtran  Devices  type  EGA-125)  mounted  on  the  tip  of  the  rod  next  to  the  model. 
Outputs  of  the  accelerometer  and  force  platform  were  monitored  by  a  dual  channel 
strip  chart  recorder  (Gould  2200).  The  net  force  on  the  rod  was  measured 
independently  and  subtracted  from  measured  forces  for  accelerated  models.  Since 
these  intertidal  animals  are  normally  attached  to  the  substratum  we  chose  to 
accelerate  our  models  parallel  to  a  boundary  1  -  2  mm  from  the  base  of  each 


model. 

Since  the  velocity  of  the  model  was  small  during  the  initial  stages  of 


acceleration  (less  than  0.1  s  after  starting  from  rest)  the  drag  on  the  body  was 
negligible  and  the  total  force  was  dominated  by  the  inertia  of  the  body  and  its 
acceleration  reaction.  Hence,  the  total  force  acting  on  the  models  can  be 
expressed  in  a  form  similar  to  eq.  4. 


=  [ft*.  V  dU/dt]  +  m  dU/dt  (12) 

where  F^f  is  the  total  force,  e*  is  the  added  mass  coefficient,  V  is  the  volume 
of  fluid  displaced  by  the  model,  m  is  the  mass  of  the  model  and  dU/djt  is  the 
acceleration.  With  measured  or  assumed  values  for  F^f,  p  ,  m,  V,  and  dU/dt  the 
added  mass  coefficient  was  determined  from  the  above  equation. 

For  all  equations  that  follow  an  inertia  coefficient  is  defined  for 
stationary  oganisms  in  accelerating  flows  as  Cm  =  1  +t<  to  include  the  effects 
of  virtual  buoyancy  (Batchelor  1967,  pg*  409),  and  values  are  shown  in  Table  1. 

Tenacity  and  Breaking  Strength  Measurements 

The  adhesive  tenacity,  <J~max,  expressed  as  the  force  per  area  required  to 
dislodge  the  organism,  was  measured  for  ei^it  species  of  intertidal  animals.  Tenacity 
measurements  for  six  of  these  species  (Collisella  digitalis,  B.  glandula,  S. 
cariosus,  Balanus  nubilus,  Mytilus  calif omianus,  and  Strongylocentrotus 
purpuratus)  were  made  on  the  exposed  coast  of  Washington  at  Tatoosh  Island  or 
Shi-Shi  Beach,  and  the  other  species  (C^  pelta,  and  Hotoacmaea  scutum)  were 
measured  at  the  Friday  Harbor  Marine  laboratory  in  Puget  Sound.  All  animals 
tested  were  solitary,  i.e.  they  were  fully  exposed  to  the  prevailing  flow  and 
were  not  shielded  by  their  neighbors.  The  outer  coast  measurments  were  conducted 
as  follows:  A  loop  of  string,  was  placed  around  the  organism  near  its  base,  or, 


for  mussels,  where  the  byssus  entered  the  shell,  and  the  string  was  pulled 
parallel  to  the  substratum  until  the  animal  was  dislodged.  The  force  of 
dislodgment  was  measured  with  a  recording  spring  scale  of  a  design  similar  in 
principle  to  that  of  Jones  and  Demetropoulos  (1968).  The  basal  area  of  each 
organism,  with  the  exception  of  M;_  calif ornianu3,  was  estimated  as  the  length- 
width  product.  The  recorded  force  was  divided  by  this  estimated  area  and 
multiplied  by  1.42  to  give  a  value  for  breaking  strength  in  shear,  CT^max,  for 
each  organism.  Multiplication  by  the  factor  of  1.42  is  necessary  because  the 
shear  stress  varies  across  the  base  of  an  organism  when  loaded  as  described 
above,  reaching  a  maximum  value  at  the  center  of  the  base.  The  maximum  stress, 
Ihe  stress  at  which  the  material  would  break,  is  1.33  or  1.50  times  the  average 
stress  (force/basal  area)  depending  on  whether  the  base  is  circular  or 
rectangular  (Timoshenko  and  Gere  1972).  Since  the  animals  measured  here  are 
generally  intermediate  between  these  two  shapes  an  average  value  of  1.42  was 
chosen.  The  shapes  of  the  organisms  measured  in  this  manner  make  this  method 
possible  by  preventing  the  string  loop  from  slipping  off  when  force  is  applied 
(eg  the  overhanging  apex  of  the  shell  of  digitalis).  The  two  remaining 
limpet  species  were  less  accommodating,  requiring  a  different  technique. 
Measurements  made  on  pelta  and  scutum  at  Friday  Harbor  were  made  using  a 
strain  gauge  force  transducer  similar  to  that  described  by  Denny  (1982).  A  short 
piece  of  wire  bent  into  a  hook  was  glued  to  the  shell  apex  using  quick-setting 
epoxy  gLue.  The  force  transducer  was  then  attached  to  the  hook  and  force  applied 
perpendicular  to  the  substratum  until  the  animal  was  dislodged.  The  force  of 
dislodgement  was  divided  by  foot  area  to  give  a  value  for  the  breaking  strength 
in  tension,  <T n>max- 

All  animals,  whether  on  the  open  coast  or  at  Friday  Harbor,  were  tested 


while  attached  to  the  natural  substratum,  and  all  were  tapped  or  otherwise 
bothered  sligvtly  immediately  before  testing  to  induce  the  animal  to  adhere  as 
best  it  could. 

The  values  of  tenacity  measured  in  this  study,  as  well  as  values  from  the 
literature,  are  summarized  in  Table  1.  With  the  exception  of  the  study  of  Miller 
(1974),  all  limpet  tenacity  values  were  determined  for  stationary  animals. 
Unfortunately  Miller's  values  for  the  moving  tenacity  of  limpets  are  somewhat 
suspect  since  the  values  reported  for  stationary  limpet  tenacity  measured  under 
similar  conditions  in  the  laboratory  (i.e.  for  animals  adhering  to  acrylic 
plastic)  are  substantially  lower  than  the  values  of  other  authors  determined  in 
the  field.  Assuming  that  the  ratio  of  moving  to  stationary  tenacitites  as 
measured  by  Miller  (1974)  are  accurate,  we  have  estimated  the  tenacity  of  moving 
limpets  by  multiplying  the  field-determined  stationary  tenacity  values  for  each 
species  by  the  average  ratio  of  moving  to  stationary  tenacity  for  all  species 
using  ditaxic  waves,  0.334  (as  determined  by  Miller  1974). 

The  breaking  strength  of  the  hydrocoral  Millepora  complanta  was  determined 
by  three-point  bending  tests.  Fresh  corals  were  collected  from  depths  of  <  1  m 
either  at  Galeta  Point,  or  from  the  algal  ridge  off  of  the  Limon  island  group, 
San  Bias  Islands,  Republic  of  Panama.  The  corals  were  cut  into  strips 
approximately  1  cm  wide  vising  a  diamond  saw;  strips  were  cut  parallel  to  the 
long  axis  of  each  coral  blade.  For  each  test  the  sample  was  simply  supported 
near  its  ends,  the  distance,  D,  between  supports  measured,  and  a  bucket 
suspended  by  a  wire  from  the  centerline  of  the  sample.  Water  was  added  to  the 
bucket  and  the  force  acting  on  the  sample  thus  increased  until  the  sample  broke. 
The  cross-sectional  dimensions  of  the  sample  were  measured  at  the  point  of 
failure  and  the  modulus  of  rupture  (a  measure  of  strength  in  bending)  calculated 


according  to  Wainwright  et.  al.  (1976), 

modulus  of  rupture  =  D  F  y/  41  (1 3) 

where  F  is  the  mass  of  the  bucket  and  water  multiplied  by  the  acceleration 
of  gravity  (9.81  m/s^),  y  is  1/2  the  thickness  of  the  sample,  and  I  is  the 
second  moment  of  area  of  the  cross-section  (=  [2/3]  width  yp)  for  a 
rectangular  cross-section.  The  modulus  of  rupture  for  the  39  samples  tested  was 
1 .23  - 1 0?  N/m^  +  1.47*10^  N/m^  (95$  confidence  intervals). 

Dimensions  of  Mussels 

Mussels  were  collected  for  morphological  measurements  from  a  0.25  m  square 
quadrat  haphazardly  thrown  onto  the  mussel  bed  at  Tatoosh  Island.  The  length  of 
each  of  these  mussels  was  measured  to  the  nearest  0.2  mm  with  vernier  calipers, 
and  their  volume  was  measured  to  the  nearest  0.1  ml  by  immersing  each  mussel  in 
a  graduated  cylinder  containing  water. 

Fifty  mussels  of  various  sizes  were  pulled  off  the  exposed  rock  surface  at 
Tatoosh  Island,  the  force  of  dislodgement  being  measured  as  described  above.  A 
plot  of  dislodgement  force  vs.  length  of  the  mussel  (Figure  3)  shows  that  the 
adhesive  tenacity  increases  in  proportion  to  l'  *  •  A  length  vs.  volume  plot  for 
182  mussels  from  Tatoosh  (Figure  4)  shows  that  volume  increases  as 
Assuming  that  projected  area  is  proportional  to  the  2/3  power  of  volume,  area 
increases  as  Thus,  at  least  to  the  accuracy  of  this  estimate, 

adhesive  tenacity  increases  at  the  same  rate  as  the  projected  area,  as  assumed 
for  the  model  discussed  later. 

If  a  mussel  is  treated  as  a  streamlined  body  with  a  thickness  of 


approximately  half  its  length,  its  about  0.2  when  end  on  to  the  flow  and 

0.8  when  broadside  at  the  Re  expected  during  wave  flows  (Hoemer,  1965),  and 
these  values  are  used  here.  is  approximately  2.0  when  the  animal  is  oriented 
broadside  to  the  flow  and  1.2  when  end-on  (Daniel  1982).  The  end-on,  minimum 
value  of  projected  area,  is  used  to  calculate  an  "equivalent  stress", 

(force  of  dislodgement/Ap>mjn)  (Table  1). 


WATER  FLOW  IN  WAVE-SWEPT  HABITATS 

Before  discussing  the  theory  of  size  limitation  in  wave-swept  organisms  it 

» 

is  useful  to  examine'  in  general  terms  the  nature  of  the  velocities  and 
accelerations  caused  by  waves. 

Intertidal  organisms  are  exposed  to  breaking  waves  and  post-breaking  flows. 
The  precise  water  movements  occurring  during  the  breaking  of  a  wave  are 
exceedingly  complex,  and  mathematical  analyses  are  only  in  the  early  stages  of 
development  (e.g.  Cokelet  1977,  1979,  Peregrine  1983,  Longuet-Higgins  1982, 
Stive  1980).  It  is  possible  however  to  draw  a  general,  approximate,  picture  of 
■these  flows.  A  wave  approaching  shore  breaks  when  the  height  of  the  wave  above 
still  water  level  is  approximately  equal  to  the  depth  of  the  water  (Galvin 
1972).  The  wave  can  break  in  a  number  of  ways,  depending,  in  part,  on  the  slope 
of  the  bottom  leading  to  the  shore.  With  gradually  sloping  bottoms  waves  curl 
over  and  the  crest  plunges  forward  at  breaking;  these  conditions  apply  for  the 
reef  at  Galeta  Point.  The  bottom  off  Tatoosh  Island,  where  much  of  this  study 
was  conducted,  slopes  steeply  down  to  deep  water  and,  consequently,  the  waves 
were  generally  of  the  type  termed  "spilling',  by  Galvin  (1 972).  ?or  this  sort  of 
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breaking,  the  wave  crest  spills  down  the  front  of  the  wave  body  rather  than 
arching  over  ahead  of  the  body  of  the  wave.  For  either  case,  at  the  moment  of 
breaking  the  wave  crest  travels  with  a  velocity,  U,  parallel  to  the  bottom, 
given  as  a  first  approximation  by  linear  wave  theory  as: 

U  =  [g  (H  +  Z)]1/2  (14) 

where  g  is  the  acceleration  due  to  gravity  (9.8  m/s^)  and  H  and  Z_  are  the  wave 
height  and  water  depth  respectively  (Weigel  1964).  Thus,  for  a  five  meter  wave 
breaking  in  water  five  meters  deep  the  crest  velocity  is  approximately  10  m/s. 

In  areas  such  as  Tatoosh  Island  where  the  bottom  slopes  steeply  away  from  the 
shore,  waves  often  break  very  near  or  directly  onto  the  emersed  rocks  of  the 
shore.  In  such  a  case,  water  with  velocity  into  the  rock  is  deflected  by  the 
impermeable  surface,  and,  except  for  the  area  immediately  under  the  wave  where  a 
plunging  crest  might  first  hit,  the  primary  flow  is  a  surge  up  the  rock  face  (a 
turbulent  bore,  or  subsequent  run-up)  at  a  velocity  initially  roughly  equal  to 
that  of  the  crest  velocity  as  calculated  above  (Carstens  1968,  Peregrine  1972). 
Water  velocity  in  post-breaking  flows  may  increase  somewhat  as  a  turbulent  bore 
collapses  at  the  shoreline  (Keller  et.  al.  I960). 

We  envision  three  different  flow  regimes  to  which  intertidal  organisms  may 
be  exposed.  They  may  experience  (1 )  periodic  flows  under  shallow  water  waves, 

(2)  a  jet  of  water  normal  to  the  rock  wall  in  the  immediate  area  of  a  plunging 
crest,  and  (3)  an  impulsive  onrush  of  a  wave  front  moving  parallel  to  the 
surface  of  the  rock  wall.  Carstens  (1 968)  described  the  flow  forces  associated 
with  the  first  two  cases.  He  found  that  under  shallow  water  waves, 
accelerations  are  relatively  small,  on  the  order  of  1  m/s^,  and  that  flow  forces 
depend  largely  on  the  fluid  velocity.  It  is  important  to  note,  however,  that 


this  estimate  of  acceleration  applies  only  to  flows  under  non-breaking  waves  and 
should  not  be  used  as  an  estimate  of  fluid  accelerations  once  a  wave  has  broken 
and  run  shoreward.  Carstens  (1968)  also  explores  the  pressures  associated  with 
the  impingement  of  a  jet  on  an  impermeable  surface.  These  impact  pressures  may 
be  quite  large,  with  values  exceeding  1.5  x  10^  Pa.  Organisms  exposed  to  such 
pressures  may  fail  by  implosive  rupture  rather  than  being  sheared  off  the  rock 
wall.  The  magiitude  of  impact  pressures  is  a  function  of  the  acoustic  velocity 
of  water  which,  in  turn,  depends  quite  strongly  on  the  fractional  volume  of 
entrained  air.  With  10#  air  by  volume,  the  acoustic  velocity  decays  to  about 
10#  of  its  value  with  no  entrained  air  (Carstens,  1968).  Since  Longuet-Higgins 
and  Turner  (1974)  estimate  an  air  content  of  about  10#  for  spilling  breaking 
waves,  impact  pressures  may  be  nearly  an  order  of  magnitude  lower  than  that 
reported  by  Carstens  (1968).  Our  records  of  wave  force  show  that  maximum  force 
is  generally  reached  50-100  ms  after  initial  contact  with  the  water.  These 
forces  are  clearly  not  due  to  impact  pressures. 

The  third  flow  regime  is  the  principal  topic  of  this  paper.  When  an 
unsubmerged  organism  encounters  the  onrush  of  a  water  front  moving  parallel  to 
the  surface,  fluid  accelerations  can  be  quite  large  (?ig.  5).  Because  the  flow 
is  not  directed  into  the  substratum  there  will  be  no  substantial  impact  pressure 
and  the  flow  is  likely  to  be  quite  turbulent.  Local  surface  topography,  water 
turbulence,  and  the  formation  along  the  substratum  of  a  layer  of  fluid  whose 
motion  is  retarded  (the  boundary  layer)  all  have  significant  effects  on  the  flow 
field  encountered  by  an  organism.  In  order  to  explore  the  magnitude  of  the 
forces  associated  with  this  flow  we  make  one  crucial  assumption:  the  fluid  speed 
and  acceleration  at  any  instant  in  time  do  not  change  as  a  function  of  distance 
from  the  substratum.  Thus,  we  assume  the  kinematics  of  the  flow  encountered  by 
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an  organism  to  be  independent  of  body  size.  Clearly,  there  are  situations  in 
which  this  assumption  is  violated  as  will  be  discussed  in  a  later  section  of 
this  paper.  However,  we  believe  this  assumption  to  be  a  valid  first 
approximation  for  solitary  organisms.  This  conclusion  is  based  on  the  following 
reasoning: 

1 )  The  establishment  of  a  boundary  layer  requires  time.  Is  there 
sufficient  time  during  the  onrush  of  breaking  wave  for  viscosity  to  result  in 
the  build-up  of  a  significantly  large  boundary  layer?  This  question  is  not 
relevant  to  "steady"'  turbulent  flows,  and  the  problem,  here,  is  to  estimate  the 
growth  in  time  of  a  turbulent  boundary  layer.  We  use  an  example  provided  by 
Schlichting  (1979):  the  growth  of  a  turbulent  boundary  layer  from  the  leading 
edge  of  a  flat  plate.  The  boundary  layer  thickness,  £  ,  is  defined  as  the 
distance  above  the  substratum  at  which  the  velocity  is  99$  of  free  stream 
velocity.  Schlichting  gives  an  equation  which  relates  boundary  layer  thickness 
to  the  distance  of  a  point  from  the  leading  edge  of  a  flat  plate: 

S  =  0.37  X  (U  x/v)0-2  (15) 

where  x  is  the  distance,  U  is  the  free-stream  velocity,  and  U  is  the  kinematic 
viscosity.  If  we  consider  a  bulk  of  fluid  moving  with  speed  U  towards  a  point 
on  a  flat  rock  wall,  we  may  estimate  the  growth  of  the  boundary  layer  in  time  at 
that  point  by  equating  x  to  Ut,  where  t  is  the  elapsed  time.  Thus,  eq.  15 
becomes: 

$  =  0.37  U0-6  t°-8  °-2  (16) 

Using  Figure  5  to  provide  numerical  estimates  of  the  values  in  the  above 
equation  and  assuming  a  kinematic  viscosity  of  10~^  m2/s  we  calculate  that  the 
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boundary  layer  would  build  up  to  its  maximum  thickness  of  only  5  x  10~3  m  ( t  ~ 
0.05  s,  U  -  5  m/s).  Of  course  the  time  averaged  value  of  &  encountered  by  an 
organism  would  be  lower  than  this  estimate  of  maximum  5  .  Thus  the  height  above 
the  substratum  at  which  water  velocity  reaches  99$  of  free  stream  is  only  a 
fraction  of  the  height  of  most  organisms  encountered  in  this  study. 

2)  A  velocity  gradient  exists  within  a  boundary  layer.  The  crucial  feature 
for  our  analysis  is  not  how  thick  the  boundary  layer  is  but  whether  there  is  a 
significant  change  in  the  velocity  as  an  organism  grows  taller.  While  for 
similar  free-stream  velocities  turbulent  boundary  layers  are  thicker  than 
laminar  ones,  they  have,  in  general,  blunter  profiles  -  the  velocity  gradient  is 
steepest  in  a  region  confined  to  a  small  distance  from  the  boundary.  Prom  this 
region  outwards,  the  change  in  velocity  is  quite  small.  At  only  5$  of  the  way 
up  into  the  boundary  layer,  the  local  velocity  is  already  greater  than  50$  of 
the  free  stream  velocity.  Prom  2Q$  of  the  boundary  layer  thickness  outwards, 
the  velocity  changes  only  by  a  factor  of  about  25$.  (See  Pig.  18.5  in 
Sclichting  1979).  Consequently,  our  assumption  of  a  constant  velocity  is 
reasonable  for  turbulent  boundary  layers  even  as  thick  as  1  or  2  cm.  Thus, 
owing  to  the  very  short  duration  of  the  initial  surge,  and  to  the  nature  of 
turbulent  boundary  layer  profiles  our  assumption  appears  valid  that  the  flow 
encountered  by  macroscopic  solitary  organisms  is  independent  of  body  size. 

The  force  trace  from  one  large  breaking  wave  is  shown  in  Pigure  5,  along 
with  the  water  velocities  and  accelerations  calculated  to  be  responsible  for  the 
force.  The  velocity  at  maximum  acceleration  cannot  be  precisely  predicted  for 
waves  in  general  and  indeed,  for  the  thirty  waves  for  which  calculations  were 
made,  no  significant  correlation  exists  between  maximum  acceleration  and  the 
velocity  at  which  that  acceleration  occurred  (r  =  0.022,  N  =  50,  P  >  0.5.) 
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MODELS 


Here  we  examine  the  hydrodynamic  forces  that  act  on  a  simple  general  model 
of  a  wave-swept  organism  as  shown  in  Figure  6.  The  model  "organism",  is  a 
hemisphere  of  radius  Lj.  attached  to  a  planar  surface.  The  passage  of  a  wave 
causes  water  to  flow  along  the  surface  and  past  the  organism,  placing  a  number 
of  forces  on  the  body.  In  the  direction  of  flow  (i.e.  parallel  to  the  surface 
plane)  there  is  a  force,  which  has  two  components:  Fjj,  the  drag  due  to  the 
instantaneous  velocity  of  the  water,  and  Fgt  the  acceleration  reaction  due  to 
the  water's  acceleration  .  By  rearranging  eq.  7  it  can  be  seen  that  F^  varies 
with  velocity,  U: 

2a  (17) 

Thus  for  the  hemispherical  model  where  =  0.5  tt  (Lj.2) 

*4  =  °*25  f  if  Lp2  u2  (18) 


The  accele rational  component  is 

Za  =  e  Cm  V  dU/dt  =  [2/5]  fC^  tr  Lj.3  dU/dt  (  n) 

where  V  is  the  volume  of  the  model  ([2/5]  tl  L^).  Cjh  for  a  hemisphere  was 
measured  to  be  1.35*  Both  2d  la.  are  assumed  to  be  applied  at  the  centroid  of 
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the  projected  area,  i.e.  along  the  midline  and,  in  the  case  of  a  hemisphere,  a 
distance  4Lr/(3n)  up  from  the  base. 

In  addition  to  5^  and  5^  which  act  only  along  the  line  of  flow,  the  model 
experiences  a  force  at  right  angles  to  the  flow  and  away  from  the  surface,  the 
lift,  P^,  which  varies  with  water  velocity  in  a  manner  similar  to  drag: 

P~1  =  0.5  f  C~i  Afo  =  0.5  <?  C~i  Lj.2  u2  (20) 

where  Ajj  is  the  basal  area  of  the  model  (Lr^).  Recall  that  lift  is  not  centered 
on  area  Afc>  but  rather,  acts  as  a  point  force  some  distance  downstream  of  center 
(the  center  of  pressure).  Lift  in  hemispheres  attached  to  walls  has  not  been 
sufficiently  studied  to  allow  precise  location  of  the  center  of  pressure,  but 
assuming  that  the  hemispherical  model  behaves  similarly  to  the  conical  models 
actually  measured,  the  lift  acts  through  a  point  along  the  midline  and 
approximately  0.24  Lj.  downstream  of  center.  As  shown  in  Appendix  2,  minor 
variations  in  the  exact  point  of  application  of  lift  do  not  substantially  affect 
the  calculations  made  here. 

If  the  model  is  to  remain  stationary  and  attached  to  the  surface  plane,  the 
forces  P<j,  and  ?x  must  be  resisted  by  the  adhesive  which  holds  the  model  in 
place.  The  forces  acting  on  the  adhesive  are  of  two  sorts.  P^  and  P^  tend  to 
push  the  model  along  the  plane,  shearing  the  adhesive.  The  average  shear  stress }  or  , 
is  thus  (Pg^  +  P^)/Afr.  Por  a  circular  cross-section  shear  stress  varies  across 
the  basal  area  and  reaches  a  maximum  at  the  center  (Timoshenko  and  Gere  1972): 

C5s1max  =  4  °^/3.  (21 ) 

Expanding  eq.  21: 
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s, max  ~  (4/3)  (?d  +  Za'/Ab 

=  (1  /3)p  £d  U2  +  (8/9)e  C £  Lj.  dU/dt 

=  Ki  U2  +  dU/dt  (22 

The  size  of  the  model,  expressed  in  this  case  as  its  radius,  L rt  thus 
affects  the  maximum  shear  stress  by  affecting  the  acceleration  reaction;  for  a 
given  velocity  and  acceleration,  the  larger  the  model  the  greater  the  shear 
stress  on  its  basal  adhesive. 

In  addition  to  a  tendency  to  slide  sideways  in  response  to  a  shear  force, 
the  model,  as  shown  in  ?igure  6,  tends  to  be  overturned,  usually  in  a  direction 
such  that  its  upstream  edge  is  lifted  and  its  downstream  edge  is  forced  down. 
This  tendency  is  due  to  the  moment  applied  by  the  shear  force  acting  some 
distance  above  the  base.  This  overturning  moment  is  counteracted  to  some  degree 
by  the  moment.  Mg,  due  to  the  lift  force  acting  downstream  of  the  center  of  the 
model.  Unless  M-|  =  there  is  some  net  tendency  to  rotate  the  model  about  the 
neutral  axis.  The  neutral  axis  is  a  line  in  the  plane  of  the  base  perpendicular 
to  the  direction  of  flow  and  passing  through  the  center  of  the  base.  Depending 
on  its  direction,  this  rotation  compresses  or  extends  the  adhesive  on  the 
upstream  half  of  the  base  and  does  the  opposite  to  the  downstream  half.  The 
adhesive  at  the  neutral  axis  is  neither  compressed  nor  extended.  The  net  moment 
(M-j  -  M2)  thus  exerts  a  stress  on  the  adhesive  in  a  direction  normal  to  the 
basal  plane.  The  magnitude  of  this  normal  stress,  Orn,  varies  with  the  distance, 
d,  from  the  neutral  axis  (Timoshenko  and  Gere  1972): 


CTn  =  (M!  -  Mg)d/I 
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v.here  I_  is  the  second  moment  of  area  about  the  neutral  axis.  ?or  a  circular 
cross-section  as  in  this  model  I_  =  0.25  tt  LjA  The  normal  stress  thus  reaches  a 
maximum  at  the  leading  and  trailing  edges  of  the  model  where  <1  =  Lj..  In  addition 
to  the  normal  stress  due  to  the  net  moment  there  is  a  normal  stress  placed 
directly  on  the  adhesive  by  the  lift.  This  normal  stress  tends  to  change  the 
location  of  the  neutral  axis  within  the  model  to  a  position  downstream  of  center 
(Timoshenko  and  Gere  1972)  but  unless  is  large  compared  to  (which  it  is 
not  for  the  bodies  we  studied)  the  effect  is  minor  and  has  not  been  included  in 
this  analysis.  Thus  the  maximum  normal  stress,  is  : 

O'n.max  -  +  (ELl_  *  Mlr/l  <24> 

Expanding  the  terms  in  eq.  2  and  rearranging: 

tf^max  =  aSfCi  U2  -  -472  f  C*  U2  +  (4/3t»)ec^  U2 
+  (32/9tt  )e  Cg,  Lj.  dU/dt 

=  K5  e  U2  -  ]i4eU2  +  K5  eU2  +  *6  dU/dt 
=  -  K4  +  K5)  +  K6fLr  dU/dl  (25) 

As  with  the  maximum  shear  stress,  the  maximum  normal  stress  for  a  given 
velocity  and  acceleration  is  affected  by  the  size  of  the  model;  the  larger  the 
radius,  the  greater  Cfc>max.  rearranging  equations  22  and  25  it  can  be  seen 
that  L-  max,  a  measure  of  the  maximum  size  of  the  model,  is  for  failure  due  to 


shear  stress: 


hr, max  =  (  ^3,  max  ~  M2)/(K 2  e  dU/dt)  (26) 

or  for  failure  due  to  normal  stress: 

jtr.max  =  (^n.max  -  U2  e  (K3  -  +  ]§))/(&e  dU/dt)  (27) 

In  the  second  case  the  adhesive  will  in  all  likelihood  fail  in  tension,  most 
materials  being  much  stronger  in  compression  than  in  tension  (Wainw right  et.  al. 
1976).  The  maximum  size  of  the  model  (or  an  organism  similar  to  the  model)  thus 
depends  on: 

1 )  The  strength  of  the  structure  or  the  adhesive  holding  the  structure  in 
place, 

2)  the  velocity  of  the  water,  and 

3)  the  water's  acceleration. 

These  parameters  may  be  empirically  measured  and  used  to  examine  critically 
the  maximum  sizes  of  organisms  inhabiting  wave-swept  environments. 

The  various  constants  derived  here  for  a  hemisphere  (K^  to  Kg)  may  easily 
be  derived  for  other  shapes  such  as  cones  and  cylinders,  and  these  values  are 
shown  in  Table  2. 

For  certain  organisms,  some  of  the  forces  calculated  above  are  negligible 
and  for  these  organisms  a  simplified  model  suffices.  For  example,  consider  a 
plate-like  organism  such  as  the  fire  coral  Millepora  complanata,  extending 
perpendicularly  out  of  a  planar  surface  •  If  the  height  of  the  organism,  L^,  is 
much  greater  than  its  depth,  the  length  of  the  moment  arm  (M  =  L^/2)  is  such  that 


the  normal  forces  exerted  at  the  edges  of  the  base  far  exceed  the  shear  force. 

The  shear  force  may  thus  be  neglected.  Similarly,  for  such  an  organism,  Cj  is 
approximately  zero  (Hoerner  and  Borst  1975).  Further  the  organism  has  a  very  small 
area  over  which  lift  can  act.  Thus,  we  neglect  lift  forces  in  the  calculation 
of  size: 


=  0. 5  C^j  W  L ^2  p  Qjjj  W  T_  L^3  dU/ dt 


(28) 
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where  W  and  T  are  the  ratio  of  the  width  of  the  plate  and  thickness  of  the 
plate,  respectively,  to  L^.  Continuing  the  calculation: 

Q'n.max  °  Bz/I 

=  3  (p/2T2)  U2  +  3  (P/T)  Ijj  dO/dt  (29) 

where  to  is  the  bending  moment,  is  the  distance  from  the  neutral  axis  (1  /2  the 
thickness)  and  1^  is  the  second  moment  of  area  ([1/12]  WT5  Lh  4).  Similar 
calculations  may  be  made  for  organisms  shaped  like  long,  thin  cylinders  and  the 
appropriate  constants  are  summarized  in  Table  2. 

One  final  and  important  addition  must  be  made  to  these  models  before  they 
can  be  productively  applied  +o  "real  world",  situations.  As  stated  previously, 
the  maximum  size  of  an  organism  in  a  certain  flow  situation  is  a  function  of  the 
strength  of  some  critical  structure  in  the  organism.  If  all  individuals  of  a  . 
oertain  species  have  the  same  breaking  strength,  all  have  the  same  mechanical 
limitation.  However,  in  any  actual  species  some  distribution  of  breaking 
strengths  is  to  be  expected,  and  for  the  organisms  examined  here,  this 
distribution  was  found  to  be  indistinguishable  from  a  normal  distribution 
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(graphical  method,  Zar,  1974).  Thus,  the  breaking  strength  distribution  of  each 
species  may  be  characterized  by  a  mean  and  a  standard  deviation,  and  for  a 
population  of  organisms  the  fraction  having  a  breaking  strength,  cr ,  of  at  least 
a  certain  magnitude  can  be  specified: 


Pg  =  1  -  fs"1  (2tt)-1/2  exp-[(a--or)2/2  s2]  d  cr 


where  s  is  the  standard  deviation,  cr  is  the  mean  breaking  stress  and  cr'  is  the 

breo.k'*'} 

observed^stress  (Zar,  1974).  Thus  Pg  is  the  fraction  of  the  population  that,  on 
average,  i£  not  dislodged  when  a  stress,  cr',  is  applied  to  each  organism. 
Conversely,  the  fraction  that  does  not  have  a  certain  adhesive  strength,  P^,  (= 

1  -  Pg)  is  a  measure  of  the  fraction  of  the  population  which  is  dislodged  try  one 
particular  wave  with  a  given  set  of  flow  parameters. 

Alternatively,  Pg  and  P^  may  be  viewed  as  an  estimate  of  the  probability  of 
survival  or  dislodgement  respectively  for  an  individual  chosen  at  random  from 
the  population.  Because  natural  selection  acts  on  the  individual  rather  than  a 
population,  this  second  viewpoint  is  more  relevant  to  the  argument  at  hand  and 


is  maintained  throu^iout  this  discussion. 

The  cumulative  probability  of  survival  or  dislodgement,  P, 


.£s,cum  or  ^d,cum’ 


respectively,  can  be  calculated: 


£3, cum  =  (^3,1  x  -s,2  x  ‘  ‘  ‘  ^s,N^ 


p 

•id,  cum 


~  ^  ”  (Zs,cum^ 


where  Pg^  is  the  probability  of  not  being  dislodged  by  each  individual  wave  and 
N  is  the  number  of  waves  encountered.  Note  that  eq.  31  requires  that  the  flow 
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parameters  of  each  of  the  N  waves  he  known.  Further,  because  each  wave 
dislodges  individuals,  the  parameters  of  the  breaking  stress  distribution  for  the 
remaining  individuals  may  vary  from  one  wave  to  the  next.  Thus,  while  correct  as 
written,  the  application  of  eq.  31  is  difficult.  A  method  for  approximating  eq. 
31  is  discussed  later. 

For  an  occurrence  such  as  dislodgeraent  or  catastrophic  breakage  it  is 
unlikely  that  (without  actually  being  dislodged)  an  organism  can  ".judge"'  its 
instantaneous  position  on  its  species'  breaking  strength  distribution.  Once 
dislodged,  it  is  too  late  to  act  on  the  knowledge.  Thus  it  would  be 
advantageous,  in  an  evolutionary  sense,  for  an  organism  to  grow  only  to  such  a 
size  where,  for  given  flow  parameters,  its  cumulative  probability  of 
dislodgement  is  acceptably  small.  The  definition  of  what  constitutes  an 
acceptably  low  probability  is  dealt  with  in  the  Discussion. 

Note  that  the  probability  of  dislodgement  is  the  sum  of  two  component 
probabilities.  From  an  examination  of  equation  5  it  can  be  seen  that  in  theory 
even  in  steady  flow  (dU/dt  =  0)  a  stress  is  present  due  to  the  drag  force.  This 
size-  and  acceleration-independent  stress  can  be  used  in  eq.  30  to  calculate  a 
size-independent  base  probability  of  dislodgement,  P^.  For  an  organism  in 
accelerating  flow  there  is  a  stress  due  to  the  acceleration  reaction  that  can 
also  be  used  to  calculate  a  probability,  in  this  case,  the  accelerations! 
component  of  the  probability  of  dislodgement.  Thus 


+  la 


(32) 


For  an  organism  of  a  certain  shape,  exposed  to  a  given  set  of  flow 
parameters  and  having  been  chosen  at  random  from  a  population  with  a  certain 

r. 
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distribution  of  breaking  strengths,  it  is  possible  to  calculate  the  increase  in 
risk  concommittant  with  an  increase  in  size,  dP^/dL. 

dP^/dL  =  (dP^/doO  *  (d  cr/dL).  (33) 


For  shear  failure  (eq.  22) 

dcr/dL  =  dlKp^Lj.  dU/dt  +  fU2)/^  =  K^dU/dt 

or  for  failure  due  to  normal  stress  (eq.  25) 


d  cr/dL  =  d[U2f  (K3  -  K4  +  K5)  +  KgeL^  dU/dt  ]/dl  =  K6e  dU/dt 
assuming  dCfl/dRe  =  0.  From  eq.  30 

fa 

dP^/dcr  =  d/s-1  (2  Tl  )-1  /2  exp-[cr-8')2/2s2]  d<*/  dcr 
/  o 

=  s-^  (2  fT)-1  /2  exp~(cr-  cr)2/2s2) 


(34) 


(35) 


(36) 


Therefore,  for  shear 


dP^/dL  =  s"1  (2tf)-1/2  eXp-[(C7-a)2/2s2]  p  dU/dt ) 


(37) 


or  for  failure  due  to  normal  stresses 


dP^/di  =  s"1  (2tT  )-1  /2  exp-[(cr-cr)2/2g2]  (Kg e dU/dt) 


(38) 


Equations  37  and  38  specify  the  change  in  Pj  per  change  in  absolute  length.  Of 


more  interest  here  is  the  change  in  per  change  in  relative  length.  This  is: 

dPa/(dL/L)  =  L-dP^/dL  (39) 

For  failures  due  to  shear  stresses 

LdP^/dl  =  s_1  (2  tr  )-1  /2  exp-[((T-&)^/2s^]  (dU/dt)  (40) 

or  for  failure  due  to  normal  stresses 

=  s-1  (2Tf)”l/2  exp-[(o'-&)2/2s2]  IiCgpdU/dt  (41) 

%  examining  eqs.  39,  40  and  41  it  can  be  seen  that  L  dP^/dl,  the  size- 
specific  increment  in  risk,  is  greatest  when: 

1)  K2  or  Kg  are  large  (i.e.  when  9m*  the  inertia  coefficient,  is  large), 

2)  £  is  small  (i.e.  when  there  is  a  narrow  range  of  breaking  strength  in  the 
population), 

3)  L  is  large  (i.e.  when  the  absolute  size  of  the  organism  is  large), 

4)  and  when  the  stress  on  the  organism  is  equal  to  the  mean  breaking 
stress.  For  a  given  set  of  flow  parameters,  &=&  for  one  certain  length  and  the 
factor  exp-[(cr-cr)2/2s^]  decreases  at  lengths  either  longer  or  shorter.  If  for 
a  given  flow  regime  this  size-specific  increment  in  risk  is  large,  an  animal 
would,  by  increasing  its  size,  be  exposing  itself  to  a  substantial  increase  in 
?d  for  each  wave  encountered.  In  such  a  case  it  can  be  reasonably  hypothesized 
that  the  size  of  a  certain  species  may,  through  the  operation  of  natural 
selection,  have  been  limited  by  mechanical  factors.  The  strength  of  the 
selection  pressure  setting  a  limit  to  size  depends  on  the  relative  disadvantage 


of  increasing  as  compared  to  the  increase  in  fitness  concommittant  with  an 
increase  in  size.  A  detailed  examination  of  this  interplay  is  postponed  to  the 
Discussion.  However,  at  this  point  it  will  be  useful  to  choose  a  cutoff  value 
for  LdP^/dL  above  which  LdP^/dl  is  considered  "large"*  Eor  reasons  treated  more 
fully  in  the  Discussion  a  value  of  LdP^/ai  >  0.1  is  considered  to  be  large 
enough  to  have  limited  the  size  of  an  organism. 


APPLICATIONS  0?  THE  MODEL 

We  now  proceed  to  examine  specific  examples  of  organisms.  Eirst  we 
consider  four  cases  in  which  the  basic  ideas  of  the  model  are  born  out:  sea 
urchins,  limpets,  mussels,  and  one  species  of  hydrocoral. 

Sea  Urchins 

The  body  of  St rongylocent  rotus  purpuratus,  the  common  purple  urchin  of  the 
Pacific  Northwest,  consists  of  a  roughly  hemispherical  rigid  test  (up  to 
approximately  0.07  m  diameter)  from  which  numerous  spines  protrude.  On  fully 
exposed  shores  individuals  generally  restrict  themselves  to  shallow  "burrows" 
worn  into  the  rock  surface  and  thus,  presumably,  avoid  the  full  brunt  of  wave 
forces.  However,  on  shores  that  are  subject  to  lesser  wave  action  (i.e.  Shi  Shi 
Beach)  urchins  are  found  fully  exposed  to  the  forces  of  the  prevailing  waves, 
and  are  thus  suitable  for  the  analysis  describe  -  uere. 

Calculations  have  been  made  here  for  S.  purp.'atns  from  Shi-Shi  beach  using 
the  C(j  and  values  shown  in  Table  1,  and  the  equation  for  shear  stress  (Eq. 


Figure  7  shows  a  plot  of  vs.  acceleration  for  four  water  velocities  (1,5, 
10,  and  20  m/s).  At  10  and  20  m/s  the  urchins  are  exposed  to  high  base 
probabilities  of  dislodgement  indicating  that  at  these  velocities  a  fully 
exposed  organism  would  be  dislodged  regardless  of  its  size.  At  the  lower 
velocities  presumably  found  in  their  semi-exposed  habitats  (1-5  m/s)  the  base 
probability  of  dislodgment  is  low  (<  0.1)  but  the  size-specific  increment  in 
is  large  (>  0.3)  (Fig.  7b, c).  Thus,  if  these  animals  experience  flows  in  the 
range  of  1-5  m/s  their  size  may  be  limited  by  mechanical  constraints. 

Limpets 

Limpets  are  small  herbivorous  gastropods  with  characteristics  conical 
shells  found  attached  to  intertidal  rocks  on  wave-swept  shores.  During  low 
tides  and  when  disturbed  limpets  clamp  their  shell  tightly  against  the  rock.  At 
high  tide,  however,  when  water  flow  is  greatest,  these  animals  move  around  to 
browse  algae.  The  C^,  and  values  used  in  the  calculations  made  here, 
along  with  the  values  defining  the  distribution  of  adhesive  tenacities  are  shown 
in  Table  1 .  The  equation  for  shear  stress  (eq.  22)  was  used  for  digitalis; 
the  equation  for  normal  stress  (eq.  25)  for  pelta  and  scutum.  Figures 
8ai,bi,ci  show  the  vs.  acceleration  curves  for  the  three  limpets  used  in  this 
study  at  four  water  velocities.  It  can  readily  be  seen  that  stationary  limpets, 
firmly  clamped  down  (dashed  lines),  run  little  risk  of  being  dislodged  by  wave 
forces.  However,  when  the  calculated  breaking  stress  values  for  moving  limpets 
are  used  in  the  calculations  (solid  lines),  both  the  base  probability  of 
dislodgement  and  the  size  specific  increment  in  risk  are  predicted  to  be  much 
higher.  This  is  most  obvious  for  C.  pelta  and  N.  scutum  (see  Fig  8a  and  b). 

The  tenacity  of  0_  digitalis  is  greater  than  that  of  either  (L_  pelta  or 
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scutum,  and  its  size  i3  smaller,  hence  for  this  species  LdP^/ai  and  are  both 
lower.  Note  that  at  a  very  high  water  velocity  (20  m/s)  the  base  probability  of 
dislodgement  is  high  (>  0.35)  for  C^_  pelta  and  scutum.  At  such  high  flow 
rates  these  two  species  would  of  necessity  have  to  stop  moving  (and  thereby  stop 
feeding)  if  they  were  to  have  any  reasonable  chance  of  remaining  attached.  In 
contrast  0.  digitalis,  a  more  common  limpet  on  very  exposed  shores  at  Tatoosh 
Island,  has  a  considerably  lower  base  risk  even  at  the  most  extreme  flow 
conditions.  For  all  three  species  an  increase  in  length  above  that  seen  in 
nature  results  in  an  increase  in  L  dP^/d L  (Figures  8aii,bii,cii).  From  these 
data  we  conjecture  that  at  high  flow  rates  (>  10  m/s)  C.  pelta  and  N.  scutum  may 
be  mechanically  limited  to  sizes  where  the  stress  placed  on  their  adhesives  is 
small  enough  so  that  they  can  tolerate  the  decrease  in  adhesive  tenacity  that 
accompanies  feeding.  At  flow  rates  of  1-10  m/s  the  base  probabilities  are  low 
(<  0.1 )  and  LdPfl/dL  is  high  (>  0.1 )  and  we  hypothesize  that  C.  pelta  and  N. 
scutum  may  be  limited  in  size  by  mechanical  factors. 

Mussels 

The  common  mussel  of  the  Washington  coast,  Mytilus  californianus,  differs 
considerably  from  the  general  morphology  of  urchins  and  limpets.  Rather  than 
having  a  shell  closely  applied  to  the  substratum  with  a  large  adhesive  base,  the 
mussel  shell  protrudes  well  above  the  substratum  and  is  held  in  place  by  an 
array  of  byssal  threads.  Despite  this  difference  in  morphology  it  was  found 
that  the  adhesive  strength  of  californianus  scales  in  the  same  manner  relative  to 
body  volume  and  projected  area  as  required  for  the  model  described  here  (eq. 

22),  and  thus  the  same  calculations,  with  minor  rearrangements,  may  be  used.  In 
order  to  use  the  equations  previously  described,  the  dislodgment  force  for  M. 
californianus  is  expressed  as  an  "equivalent  stress",  by  dividing  force  by 


a  calculated  projected  area,  A^.  turn  varies  with  the  direction  of  flow 

relative  to  the  mussel's  long  axis,  and  is  at  a  minimum  (Ap,min>-in  =  0.126 
L^*91;  l  in  m)  when  flow  is  along  the  axis  and  a  maximum  (Ap>max  =  0.276  L^l ) 
when  flow  is  perpendicular  to  the  axis.  The  equations  appropriate  for 
calculating  Lmay  are  thus: 

f^eq,max  =  Ap  U2)/A p>min  +  (C  Cm  V  dU/dt)/Ap>mjn  (42) 

Ap/Ap>mjn  equals  1  when  the  mussel  is  end-on  to  the  flow  and  2.1 9  when  the 
mussel  is  broadside.  Using  the  relationship  shown  in  Figure  4  (V,  [in  ] 

=  0.061  L2-8^,  L  in  m) 

q'eq.nax  ~  0-5f  ik,  +  0-48f  C„  I  «/dt  (43) 

Values  of  are  shown  in  Figure  9  as  a  function  of  acceleration.  From 
Figure  19  it  can  be  seen  that  even  at  a  relatively  low  velocity  for  this  exposed 
environment  (10  m/s)  a  mussel  oriented  broadside  to  the  flow  has  a  very  hi$i 
value  of  L  dP^/dJj.  Consequently,  even  at  10  m/s  the  size  of  a  solitary  mussel 
with  this  orientation  would  be  severely  limited  by  mechanical  factors.  At 
higher  velocities  the  base  probability  of  dislodgement  approaches  unity-  Thus, 
solitary  mussels  oriented  directly  broadside  to  the  prevailing  flow  appear  to  be 
mechanically  unlikely  beasts,  and  indeed  they  are  very  seldom  found.  Mussels 
oriented  end-on  to  the  flow  present  a  more  feasible  picture  (Figure  9).  Here, 
even  at  hi$i  velocities  (20  m/s),  both  the  base  probability  and  size-specific 
increment  in  risk  are  low.  If  water  moved  only  along  one  axis  these 
appropriately  oriented  mussels  would  be  at  little  risk  from  the  flow;  however, 
in  breaking  waves  the  unpredictable,  turbulent  nature  of  flow  renders  accurate 
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orientation  unlikely  except  in  surge  channels  where  the  direction  of  flow  is 
more  predictable.  Indeed  mussels  may  be  found  oriented  in  a  direction  parallel 
to  flow  in  such  habitats  on  Tatoosh  Island.  ?or  mussels  exposed  to  flows 
neither  end  on  nor  broadside  the  force  encountered  will  be  intermediate  between 
broadside  and  end-on  values,  and  P^  and  LdP^/dL  likewise  will  be  intermediate. 

As  with  the  organisms  discussed  above,  the  size  specific  increment  in  risk  for 
M«  californianus  increases  as  the  size  of  the  animal  is  increased  above  that 
observed  in  the  field  (Pigure  19). 

Coral 

Millepora  complanata,  one  of  the  "fire  corals"-  of  the  western  Caribbean,  is 
a  platelike  hydrocoral  that  inhabits  the  surf  zone  of  reefs.  At  Galeta  Point 
the  plates  grow  perpendicular  to  the  substratum  and  are  generally  oriented 
broadside  to  the  prevailing  flow  of  breaking  waves.  VL_  complanata  typically 
grows  at  depths  of  0-4  m  and  is  thus  routinely  exposed  to  the  large  accelerations 
and  velocities  accompanying  breaking  waves. 

For  this  coral  the  thickness  of  the  plate  at  the  base  is  0.150  of  the 
height,  1^.  Appropriate  values  for  C^  (=  1.8)  and  (=6)  have  been  taken  from 
the  studies  of  flat  plates  by  Keulegan  and  Carpenter  (1958),  and  these  values 
have  been  inserted  into  eq.  29  to  calculate  stress.  Using  these  stress  values 
the  curves  of  P^  and  L  d  P^/dl  have  been  calculated  for  various  velocities  and 
accelerations  (Pigure  10). 

While  at  velocities  high  for  this  reef  environiw.. ,  (>1 0  m/s)  the  base 
probability  of  dislodgement  is  very  high  (>  0.8)  it  is  unlikely  that  waves  high 
enou^i  to  cause  such  velocities  (>5m)  occur  at  the  site  where  I^>max  was 
measured.  Por  the  maximum  wave  heights  observed  at  Galeta  Point  (2-3m,  J.  Cubit 


pers.  com.)  velocities  of  7-5  m/s  are  more  likely  (eq.  14).  At  these  velocities 
the  base  probability  is  relatively  low  and  L  dP^A?!  *s  high  (>  0.2).  On  this 
basis,  it  seems  likely  that  the  size  of  complanata  at  Galeta  Point  may  be 
limited  by  the  mechanical  consequences  of  accelerational  flow.  At  other  sites  in 
the  Caribbean  complanata  forms  plates  in  deeper  water  (2-4  m),  but  is  present 
only  as  an  encrusting  form  in  the  surf  zone  (S.  Palumbi,  pers.  com.)  This 
difference  in  morphology  could  be  occurring  in  response  to  the  mechanical 
limitations  cited  here. 

We  now  examine  six  species  that  are  apparent  exceptions  to  the  ideas 
discussed  above. 


Snails 

Data  for  three  species  of  snails  commonly  found  on  Tatoosh  Island  (Thais 
canal iculata,  Thais  emarginata,  Littorina  scutulata)  are  shown  in  Table  1 . 
Probability  calculations  have  been  made  using  the  C^  and  values  shown  in 
Table  1  and  using  eq.  22.  Compared  to  the  tenacities  of  limpets,  the  stationary 
adhesive  strength  of  these  snails  is  quite  low,  with  the  consequence  that  at 
rapid  water  velocities  (10-20  m/s)  the  size-independent  probability  of 
dislodgement  approaches  1  (see  Figure  11).  At  more  moderate  water  velocities 
(1-10  m/s)  the  stationary  adhesive  strength  of  these  animals  is  sufficient  to 
result  in  a  lower  base  probability.  However,  in  the  cases  of  T^  canaliculata  and 
T.  emarginata,  the  decrease  in  tenacity  associated  with  locomotion  causes  the 
size-independent  probability  of  dislodgement  for  moving  animals  to  approach 
unity  for  velocities  greater  than  5  m/s. 

For  all  three  species,  at  moderate  flow  conditions  (5  m/s,  <  500  m/s2)  the 
L  dP/i/di  values  based  on  stationary  tenacity  may  be  substantial.  At  main  stream 


velocities  (10-20  m/s)  or  when  the  snails  crawl  while  the  water  velocity  exceeds 
1  m/s,  the  size-specific  increment  in  risk  is  quite  low  (<  0.01)  at  all 
accelerations  due  to  the  high  base  probability  of  dislodgement.  Consequently, 
for  the  conditions  that  would  be  encountered  by  a  moving  animal  exposed  to  even 
moderate  flows,  its  risk  of  dislodgement  is  governed  primarily  by  the  size 
independent  probability  of  dislodgement.  Thus,  if  the  tenacity  values  measured 
by  Miller  (1 974)  are  correct,  these  snails  are  not  capable  of  withstanding  exposure  •»<» 
the  velocities  and  accelerations  accompanying  breaking  waves.  This  prediction 
is  indeed  born  out  since  these  animals,  while  abundant  on  moderately  exposed  (T. 
emarginata)  to  fully  exposed  (T.  canaliculate  L.  scutulata)  shores,  confine 


themselves  to  microhabitats  where  they  can  avoid  the  prevailing  flow.  During 
periods  of  rapid  water  flow  T^  canaliculata  and  T^_  emarginata  on  Tatoosh  Island 
are  typically  found  in  crevices  or  nestled  into  the  interior  spaces  of  mussel 
beds,  and  scutulata  hidden  between  acorn  barnacles  and  in  crevices  in  the 
rock  surface.  It  may  be  conjectured  then  that  the  mechanical  environment 
constrains  the  behavior  of  these  snails  and  imposes  an  indirect  limit  to  their 
body  size:  they  must  remain  small  enough  to  utilize  microhabitats  not  exposed  to 
the  full  impact  of  breaking  waves  (sensu  Kohn  1971,  Emson  and  Faller-Fritsch 
1976).  Foraging  is  thus  limited  to  times  when  the  environment  is  mechanically 
'safe ' . 

Acorn  Barnacles 

The  adhesive  strength  of  three  species  of  acorn  barnacle  (Semibalanus 
cariousus,  Balanus  glandula,  and  nubilus)  were  measured  on  Tatoosh  Island 
(Table  1).  Probability  calculations  were  made  using  the  C^,  and  values  shown 
in  Table  1  and  using  eq.  22.  The  adhesive  tenacity  of  these  organisms  is  quite 
hi$i;  consequently  the  size-independent  probability  of  dislodgment  is  low  and 


the  size  specific  increment  in  risk  even  lower  (see  Figure  12  for  cariosus). 
It  thus  seems  unlikely  that  the  size  of  these  organisms  is  limited  by  the 
mechanical  factors  considered  here.  However,  there  migfrt  be  a  mechanical  limit 
to  the  size  of  the  cirral  net,  which  must  be  extended  into  the  flow  during 
feeding,  and  thereby  of  the  whole  organism. 

DISCUSSION 

2a.  LdPfl/dL  and  Natural  Selection 

We  believe  the  theory  and  results  presented  here  to  be  a  reasonable 
description  of  the  consequences  of  one  wave  interacting  with  a  population  of 
organisms.  The  conclusions  drawn,  however,  rely  on  the  argument  that  over  the 
course  of  many  waves  large  animals  are  removed  at  a  rate  sufficiently  greater 
than  that  at  which  smaller  animals  are  removed  to  cause  a  selection  pressure 
towards  smaller  size.  Thus,  we  argue  that,  above  a  certain  size,  an  increase  in 
size  results  in  a  decreased  contribution  to  the  next  generation. 

In  order  to  justify  these  conclusions  we  must  know  something  of  the 
probability  that  an  individual  will  survive  to  reproduce.  Equation  31  presents 
a  method  for  the  exact  solution  to  the  probability  of  survival,  but,  as  noted, 
is  difficult  to  apply  in  that  it  requires  a  knowledge  of  the  velocity  and 
acceleration  accompaying  each  of  the  many  waves  occurring  in  the  time  before  an 
organism  reproduces.  Our  present  knowledge  of  the  long-term  distribution  of 
velocities  and  accelerations  occurring  in  various  microhabitats  on  exposed 
shores  is  not  adequate  to  allow  for  a  precise  application  of  eq.  31*  It  is 
possible,  however,  to  use  an  approximate  form  of  eq.  31  as  a  tool  for  comparing 


the  survivorship  of  organisms  of  various  sizes. 


We  first  examine  the  case  of  a  motile  organism.  When  an  organism  such  as  a 
limpet  or  snail  moves,  the  status  of  its  adhesion  to  the  substratum  is 
constantly  changing  -  it  may  crawl  over  a  bit  of  debris  or  an  air  bubble  and 
thus  be  less  strongly  adherent,  or  it  may  crawl  over  a  clean  rock  surface  and 
adhere  quite  well.  We  assume  here  that  as  the  organism  crawls,  its  adhesive 
tenacity  always  lies  within  the  instantaneous  distribution  shown  by  its  species 
at  that  site,  but  that  the  adhesive  tenacity  of  an  individual  at  one  time  is 
independent  of  that  at  another  time.  Thus,  for  times  short  enough  that  growth 
can  be  neglected  the  probability  of  an  individual  not  being  dislodged  by  a  wave 
of  given  prameters  is  independent  of  time  and  equal  to  Pg  (eq.  31). 

Over  a  period  of  time  N  waves  strike  an  individual.  Assume  that  the 


probability  of  encountering  a  velocity  within  a  certain  range,  Uj,  is  Pjj^  for  a 


certain  wave,  and  the  probability  of  encountering  an  acceleration  within  a 


certain  range,  A-j ,  is  ,  for  the  same  wave.  Further  assume,  as  suggested  by 

1 

the  data  cited  earlier,  that  there  is  not  corelation  between  the  velocity  and 


acceleration  encountered.  Given  these  assumptions,  the  number  of  waves  with  a 


velocity  and  acceleration  within  these  ranges  is  NPy  .  Equation  31  can  thus  be 


approximated  as 


P  =  (P  ) 
Is,  cum  'c-8' 


_  _  Mu  Pa 

x  (PJ  -a  \ 


<?>) 


Mumax  ^Amax 


(44) 


for  the  total  range  of  velocities  and  accelerations  between  U-j  and  Umax  Al 
and  A mflY.  Each  Pg  is  evaluated  per  eq.  30  at  the  stress  caused  by  the  mean 
velocity  and  acceleration  of  the  range  applicable  to  each  expression.  Thus  eq. 
31  for  the  cumulative  probability  of  survival  can  be  rewritten: 


~s>cum  -  iy  L—  £A.  iii  £s 

i*i  i  j--.  — i 

If  the  probability  distributions  of  Py  and  ^  are  known,  Pg  cum  can  be 

— i  —1  — 2 — 

evaluated.  The  distribution  of  Py  for  intertidal  habitats  has  not  been 

—  V. 

measured,  but  sufficient  information  is  known  about  waves  in  both  deep  and 
shallow  water  to  arrive  at  an  estimate.  The  height  of  waves  in  deep  water 
approximately  follows  a  Rayleigh  distribution  (U.S.  Army  1977). 

I(H>H')  =  exp  -(H7Hrma)  2  (46 

where  P(h>H')  i-3  the  probability  that  the  wave  height  is  greater  than  a  value 

H',  and  H^g  is  the  root-mean-square  wave  height.  The  orbital  velocity  of  water 
within  a  wave  is  a  function  of  the  wave  height  (see  eq.  14).  If  we  assume  that 
waves  break  near  enough  to  shore  so  that  the  velocity  encountered  by  an  organism 
on  the  shore  is  equal  to  the  crest  velocity  at  breaking  (i.e.  the  bottom  slope 
is  steep),  and  that  the  wave  breaks  when  its  height  is  equal  to  0.7  of  the  water 
depth  (Weigel,  1964),  eq.  46  can  be  rewritten  as 

I(U>U')  =  exp  -  [U'2/(1.7  g  Irms)]2  (47) 

where  P(u>t/)  is  the  probability  that  the  velocity  is  greater  than  U'.  can 

-  -i 

thus  be  approximated  as 

=  I(U>U')  “  I(U>U")  (48) 


where  u*  end  jji»  are  the  ends  of  the  range  such  that  U">U'.  In  making  these 
assumptions  we  ignore  the  increase  in  wave  height  as  deep-water  waves  shoal. 

is  unknown.  Rather  than  assuming  some  probability  function,  we  set  the 


acceleration  to  a  value  of  500  m/s^  assuming  that  at  some  point  during  wave 
breaking  this  value  will  be  encountered.  This  acceleration  is  larger  than  most 
measured  on  Tatoosh  Island,  but  as  noted  in  the  Material  and  Methods  these 
measurements  probably  underestimate  the  true  values.  Equation  45  thus  becomes. 

Zs.cum  =  exP  fij  in  is  (49) 

c-i  — i  — 

where  Pg  is  evaluated  for  the  mean  velocity  within  each  range  y  and  for  an 
acceleration  of  500  m/s^.  Given  values  for  H^g,  the  hydrodynamic  coefficients 
for  the  given  organism  (C^,  C^,  Gm),  the  population  distribution  of  breaking 
strengths,  and  the  animal's  size  an  approximation  of  the  probability  of 
surviving  the  impact  of  N  waves  can  be  calculated. 

Equation  49  is  based  on  the  assumption  that  the  adhesive  tenacity  of  the 
individual  at  one  time  is  independent  of  that  at  another  time.  While  this 
assumption  seems  valid  for  moving  limpets  or  snails,  it  does  not  seem  reasonable 
for  an  organism  that  is  permanently  cemented  in  place,  such  as  an  acorn 
barnacle.  In  this  case  it  seems  more  reasonable  to  assume  that  each  organism 
has  an  adhesive  tenacity  that  is  independent  of  time.  This  time-independence 
can  be  incorporated  into  an  approximation  of  Eq.  31  as  follows:  Consider  an 
organism  of  fixed,  but  initially  unknown,  adhesive  tenacity.  Before  this 
organism  is  subjected  to  wave  forces  it  can  only  be  said  that  its  adhesive 
tenacity  lies  within  the  distribution  of  adhesive  tenacities  shown  by  its 
species  at  the  particular  site  being  examined.  After  the  organism  survives  N 
waves  information  is  available  as  to  the  organism's  adhesive  tenacity:  Each  of 
these  waves  has  subjected  the  organism  to  a  stress,  and  the  organism  has  not 
been  dislodged.  The  oganisra  can  thus  possess  less  than  a  certain  value  of 


t 

adhesive  tenacity,  &,  only  if  every  one  of  the  N  waves  had  a  velocity  and 
acceleration  such  that  the  stress  placed  on  the  organism  was  less  than 
Assuming,  as  before,  a  value  for  acceleration,  the  stress  placed  on  the  organism 
by  each  wave  is  determined  solely  by  the  water  velocity.  We  assume  that  the 
distribution  of  velocities  is  as  described  by  eq.  47*  Thus  the  probability  that 
each  of  N  waves  has  a  velocity  less  than  U|_  is 

(50) 

and  the  probability  that  an  organism  has  less  than  a  certain  breaking  stress  is 

4  n-£(p>tf)g  (5D 

?or  example,  given  an  H^g  of  1.0  m  the  probability  that  a  given  wave  has  a 
velocity  greater  than  0.1  m/s  is  very  high  (>0.99).  Thus,  1  -  P(u>u«)  is 
nearly  zero  and  [l-P(u>U')  ^  is  smaller  still.  Having  encountered 
N  waves  and  survived  there  is  a  very  low  probability  that  the  organism  has  an 
adhesive  strength  incapable  of  withstanding  the  stress  caused  by  a  velocity  of 
0.1  m/s.  Thus  the  probability  of  surviving  a  wave  with  a  velocity  of  0.1  m/s  is 
higi.  ?or  organisms  with  a  fixed  value  «f oj  the  probability  of  surviving  the 
impact  of  one  wave  after  having  already  survived  the  impact  of  N  waves  is, 
determined  by 

Zs'  =  1“(P.d  t 1  -£(U>U')]~)  (52) 

Inserting  this  adjusted  value  of  Zb'  into  eq.  49- 

1max 

21  fij.  m  Ze'  (55) 
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where  M  is  the  total  number  of  waves  encountered. 

In  this  approximation  of  Eq.  31,  Pg'  is  evaluated  for  each  wave.  An 
intertidal  organism  may  encounter  well  over  a  million  waves  in  a  year,  and 
calculating  P's>oum  for  a  year's  time  would  be  extremely  laborious.  A  further 
approximation  can  be  made  by  re-evaluating  Pg'  on  a  daily  rather  than  a  wave-by¬ 
wave  basis: 

"I"**'*  Uu 

a,  cum  =  exp  J_  N  Z  fij .  In  (1-P.j  [1-gflwj.)]  ®  (54) 

i-\  —  k.  —  - 

where  N  is  the  number  of  waves  per  day  and  D  is  the  number  of  days. 

As  approximations  of  the  exact  solution  to  survivorship  (equations  31) 

equations  49  and  54  can  be  used  to  compare  survivorship  of  organisms  of 

different  sizes.  For  present  purposes  this  is  done  by  assuming  a  constant  value 

for  H^g  and  calculating  PgtCum  or  PngtCum  83  a  function  of  time  for  organism 

of  the  same  species  but  different  characteristic  lengths.  It  should  be  noted, 

however,  that  these  calculated  survivorship  curves  are  intended  only  as  a  means 

to  compare  different  sized  members  of  a  species,  not  as  an  accurate  calculation 

of  the  actual  survivorship  on  a  real  shore.  Several  aspects  of  the 

approximation  procedure  must  be  refined  or  modified  before  accurate  calculations 

of  "real  world"'  survivorship  can  be  made.  For  example,  1 )  the  H^g  of  the  sea 

varies  widely  from  day  to  day  and  from  site  to  site,  2)  The  approximation  of  Prr 

— £ 

by  a  Rayleigh  function  based  on  deep-water  wave  he i git  is  simply  an  educated 
guess  as  to  the  true  distribution  of  flows  in  the  surf  zone  and  will  undoubtedly 
have  to  be  modified  to  apply  directly  to  any  specific  site,  3)  the  possible 
effects  of  variation  in  acceleration  have  not  been  taken  into  account,  and  4) 


the  effects  of  the  organism's  growth  have  been  ignored. 

Figures  13  and  14  show  the  results  of  applying  these  approximations  to 
several  of  the  species  studied  here.  Due  to  their  very  high  base  probability  of 
dislodgement,  and  the  likelihood  that  their  size  is  limited  by  other  than 
mechanical  factors,  the  snails  emarginata,  T.  canaliculta  and  scutulata 
are  not  included  in  this  examination. 

Animals  are  assumed  to  encounter  5000  waves/day.  In  general  a  value  of 
Irms  =  1*0  m  has  been  used,  as  suggested  for  exposed  shores  of  the  Pacific 
Northwest  by  the  U.S.  Army  (1977).  For  those  species  occurring  in  more  protected 
areas  (S^  purpuratus,  M.  complanata)  the  value  of  H^g  has  been  reduced 
accordingly  (0.5  m  and  0.25  m  respectively).  Except  where  noted  eq.  54  (which 
assumes  fixed  adhesive  tenacity)  has  been  used. 

In  order  to  examine  whether  mechanical  factors  are  important  selective 
pressures  limiting  the  size  of  organisms  in  a  particular  population  one  must 
compare  the  increase  in  reproductive  output  that  accompanies  an  increase  in  size 
with  the  decrease  in  survivorship  that  accompanies  that  increase  in  size.  A 
proper  analysis  of  this  question  requires  knowledge  of  the  life  histories 
(especially  the  age-  and  size-specific  reproductive  values)  of  the  organisms 
being  considered  as  well  as  of  the  distribution  in  time  of  waves  of  different 
sizes  impinging  on  a  particular  site.  Because  such  information  is  not  yet 
available  for  the  animals  we  are  considering  here,  we  present  the  following 
simple  analysis  in  order  to  arrive  of  a  "rule  of  thumb"’  for  deciding  whether 
mechanical  factors  are  likely  to  be  important  in  limiting  the  size  of  wave- 
swepth  organisms  (i.e.  for  designating  L  dP^/dl  as  "’large"  or  "small"'). 

For  the  sake  of  simplicity  we  assume  that  organisms  produce  young  only  once 
per  year,  and  maintain  a  constant  size  throughout  the  year.  For  organisms  which 


grow  isometrically  (as  assumed  for  the  animals  studied  here)  a  doubling  in 
linear  dimensions  results  in  an  eightfold  increase  in  internal  volume.  We 
assume  that  this  increase  in  volume  is  used  to  produce  more  young  or  young  of  a 
greater  viability.  In  either  case  it  is  assumed  that  for  a  doubling  in  linear 
dimension  eight  times  as  many  progeny  are  contributed  to  the  next  generation. 
Given  these  two  assumptions,  a  selective  pressure  towards  increased  size  should 
exist  unless  a  doubling  in  linear  dimension  results  in  a  greater  than  eightfold 
decrease  in  survivorship  to  reproduction.  In  other  words,  if  the  survivorship 
ratio  for  an  organism  of  one  size  relative  to  another  member  of  that  species 
twice  that  size  exceeds  eight,  it  is  postulated  that  a  selection  pressure 
tending  to  limit  size  can  exist  due  to  the  mechanical  constraints  examined  here. 

This  is  clearly  a  very  simplified  approach  to  the  demographics  of 
intertidal  organisms,  however  it  is  sufficient  for  the  comparative  purposes  for 
which  it  is  employed  here. 

a.  A  Sea  Urchin  -  purpuratus:  At  the  end  of  one  year  urchins  of  the 
maximum  size  observed  at  Shi  Shi  Beach  are  41  times  more  likely  to  have  survived 
than  urchins  of  twice  this  size  (Figure  13a).  Conversely,  urchins  half  the 
maximum  observed  size  have  only  a  3*8  times  greater  probability  of  survival  than 
urchins  of  the  maximum  size.  On  the  basis  on  the  criterion  cited  above,  it 
appears  feasible  that  due  to  mechanical  factors  a  selection  pressure  exists 
limiting  the  size  of  this  urchin. 

b.  A  Mussel  -  Mytilus  calif orni anus:  At  the  end  of  one  year  a  mussel 
oriented  broadside  to  the  flow  is  1 56  times  more  likely  to  survive  at  the 
maximum  size  observed  on  Tatoosh  island  than  at  twice  that  size  (Figure  13b).  A 
mussel  half  the  observed  maximum  size  is  only  4*4  times  as  likely  to  survive  as 
at  the  maximum  observed  size.  Thus  mussels  oriented  broadside  to  its  flow  could 


be  limited  by  mechanical  factors.  However,  for  mussels  oriented  end-on  to  the 
flow  the  observed  maximum  size  is  only  5-6  times  as  likely  to  survive  as  a 
mussel  twice  this  size  (Figure  13e).  Due  to  the  unpredictable  direction  of  flow 
in  natural  environments,  any  given  mussel  is  oriented  neither  precisely  end-  on 
nor  broadside  to  most  waves.  The  appropriate  survivorship  ratio  thus  lies 
somewhere  between  156  and  5*6  and  is  likely  to  be  greater  than  8.  It  seems 
feasible  that  the  size  of  solitary  calif ornianus  is  limited  by  mechanical 

factors. 

c.  A  Coral  -  Millepora  complanata:  After  only  one  month's  exposure  to 
waves  with  Hrms  of  0.25  m  a  coral  blade  of  the  maximum  size  observed  on  the 
Galeta  reef  is  90  times  more  likely  to  remain  intact  than  a  blade  twice  the 
observed  size  (Figure  13c).  At  the  end  of  a  year  this  ratio  is  1700.  In 
contrast  to  mussels  and  urchins,  however,  a  coral  blade  that  is  half  the 
observed  maximum  size  is  approximately  16  times  as  likely  to  survive  as  is  a 
blade  of  the  maximum  observed  size  and  one  might  be  led  to  conclude  that  a 
selection  pressure  exists  for  this  coral  to  be  smaller  than  actually  observed. 

Two  factors  must  be  taken  into  account  when  examining  these  results:  1 ) 

The  reproductive  output  of  complanata  (when  reproducing  sexually)  is  probably 
more  a  function  of  the  volume  of  polyps  rather  than  of  the  volume  of  the  entire 
plate.  Assuming  that  polyps  stay  the  same  size  but  increase  in  number  as  the 
plate  grows,  the  reproductive  output  of  the  colony  should  increase  in  proportion 
to  the  surface  area  of  the  plate,  i.e.  a  fourfold  increase  for  a  doubling  in 
length.  2)  Unlike  mussels  and  urchins  which  in  all  likelihood  die  if  dislodged 
(Dayton  1973),  broken  blades  of  com  plan  ta  occasionally  land  on  favorable 

substratum,  re-attach  and  grow  (pers.  obs.).  Breakage  may  thus  serve  as  a 
dispersal  mechanism  (Highsmith  1982).  These  two  factors  lead  to  different 
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possible  interpretations  of  the  results.  The  small  increase  in  reproductive 
output  with  an  increase  in  size  for  this  colonial  or^nism  (relative  to  an 
individual  organism  such  as  an  urchin)  would  place  a  stronger  selection  pressure 
against  increased  size  if  breakage  results  in  death.  However,  if  an  increased 
probability  of  breakage  represents  in  an  increased  probability  of  dispersal, 
and,  if  dispersal  is  sufficiently  advantageous,  a  selection  pressure  could  work 
to  increase  the  mean  size  of  a  population.  Without  further  data  concerning  the 
survivorship  and  fecundity  of  broken  blades,  the  role  of  mechanical  factors  in 
limiting  the  size  of  M.  complanata  remains  questionable. 

d.  An  Acorn  Barnacle  -  cariosus:  An  individual  cariosus  of  the 
maximum  size  observed  on  Tatoosh  Island  has  only  an  approximately  6$  greater 
chance  of  surival  than  an  individual  twice  this  size  (Figure  13d).  In  this  case 
the  mechanical  factors  described  here  clearly  do  not  play  as  important  role  in 
size  limitation. 

e.  Limpets  -  0_  digitalis,  (L_  pelta,  TL_  scutum:  The  survivorship 
curves  for  the  three  limpet  species  are  shown  in  Figure  14.  The  adhesive 
tenacity  of  stationary  limpets  of  all  three  species  is  so  great  that  mechanical 
limitation  to  size  does  not  seem  feasible  (Figures  14aii,  bii,  cii).  The 
decrease  in  tenacity  accompaying  locomotion,  however,  leads  to  another  picture 
(Figures  14ai,  bi,  ci).  In  these  cases  eq.  49  (which  assumes  randomly  varying 
adhesive  tenacity)  has  been  used  to  calculate  Pg ? cum.  Because  the  predicted 
rate  of  dislodgement  of  a  population  of  moving  snails  is  high,  survivorship  is 
plotted  against  number  of  waves  rather  than  months.  For  all  three  species, 
after  only  a  few  waves  a  survivorship  ratio  of  at  least  eight  is  reached  between 
the  maximum  observed  size  and  twice  this  size.  For  example,  after  50  waves  the 
ratio  for  C.  pelta  is  92  and  for  N.  scutum,  455-  After  150  waves  it  is  13  for 
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C.  digitalis.  The  survivorship  ratios  Lj.:  2Lr  and  1^/2:  Lp  increase  with  the 
number  of  waves,  the  latter  exceeding  8  within  the  first  70  waves  for  IL_  scutum 
and  pelta  and  within  350  waves  for  digitalis. 

These  survivorship  data  serve  to  emphasize  the  conclusion  reached  earlier. 
j*br  limpets  there  is  a  large  risk  associated  with  movement.  The  smaller  the 
limpet  the  less  this  risk,  and  based  on  these  approximate  calculations,  for  at 
least  the  first  30-70  waves  during  a  period  of  movement  (150  to  350  waves  for  C. 
digitalis)  the  survivorship  ratio  is  such  that  the  mechanical  constraints 
discussed  here  are  feasible. 

The  general  picture  emerging  from  the  approximate  survivorship  curves 
presented  here  is  one  of  size  limitation  throu^i  catastrophe.  Where 
survivorship  is  high  at  the  end  of  a  year  (>  0.90)  survivorship  ratios  are 
necessarily  low  (i.e.  £!_  cariosus,  M.  californianus  end-on).  Conversely,  where 
survivorship  is  low  (<  0.10)  survivorship  ratios  exceed  the  value  of  eight 
required  to  invoke  a  selection  pressure  for  smaller  size.  Thus,  if  the 
mechanical  limits  to  size  proposed  here  are  operative  in  nature,  they  will  have 
their  most  profound  effect  under  conditions  where  within  one  reproductive  cycle 
a  large  proportion  of  a  population  is  dislodged.  A  winter  with  exceptionally 
rough  storms  would  constitute  such  a  condition  for  many  species.  Alternatively, 
such  conditions  could  apply  annually  to  an  organism  with  a  short  maturation 
time,  one  which  matures  during  summer's  calm  seas  but  may  be  ravaged  by  an  early 
fall  storm  (e.  g.  Postelsia  palmaeformis). 

Further,  the  survivorship  curves  presented  here  suggest  an  approximate 
"critical  value"'  for  Ld  P^/dl.  The  two  examples  in  which  survivorship  ratio  at 
the  end  of  a  year  is  less  than  eight  (&_  cariosus,  M.  californianus  end-on)  have 
LdPj/dl  values  less  than  0.10  at  the  maximum  size  observed  in  nature.  In  all 


cases  where  the  survivorship  ratio  exceeds  8  at  the  end  of  a  year  LdP^/dL 
exceeds  0.10.  It  is  on  this  basis  that  we  have  chosen  0.10  as  the  cut-off  value 
for  designating  LdP^/dl j  as  "large"'  or  "small".  This  value  is  intended  only  to 
be  a  "rule  of  thumb".  Any  rigorous  definition  of  a  critical  value  for  LdP^/dL 
above  which  size  is  limited  by  mechanical  factors  and  below  which  size  is 
limited  by  biological  factors  must  take  into  account  both  the  long-term  flow 
regime  at  a  particular  site  and  the  life  history  and  demographics  of  the 
organism  in  question. 

The  examples  cited  above  suggest  that  the  approach  used  in  this  study  may 
provide  a  useful  tool  for  examining  the  size  of  solitary  inflexible  organisms 
exposed  to  flows  of  high  velocity  and  high  acceleration.  In  some  cases 
mechanical  factors  alone  may  form  an  adequate  explanation  for  the  observed  size 
limits.  In  other  cases  mechanical  factors  are  apparently  of  minor  direct 
importance  in  size  limitation,  althou^i  an  indirect  affect  (as  in  the  case  of 
the  snails  discussed  above)  is  still  possible.  In  these  cases  a  search  should 
be  made  for  a  biological  limit  to  size,  or  alternatively  for  a  mechanical 
mechanism  setting  an  optimum  size  less  than  the  maximum  size  possible. 

Assumptions  of  the  Hypothesis 

While  this  approach  may  be  useful  in  many  cases,  we  do  not  intend  to  imply 
that  it  can  be  applied  to  all  wave-swept  organisms,  nor  that  it  is  a  complete 
description  of  the  forces  those  animals  encounter.  In  particular,  there  are 
several  aspects  of  the  flow  forces  associated  with  waves  for  which  the  model 
does  not  account,  and  several  of  the  model's  assumptions  require  further 
discussion. 


1.  The  derived  equations  do  not  include  a  term  for  a  moment  tending  to 


rotate  an  organism  about  an  axis  perpendicular  to  its  basal  plane,  and  thus 
torsional  forces  are  ignored.  Such  torsional  stresses  would  arise  if  5^  were  to 
act  along  a  line  lying  to  one  side  of  the  center  of  adhesive  area.  It  seems 
unlikely  that  this  happens  in  organisms  that  are  approximately  radially 
symmetrical,  such  as  limpets  and  barnacles,  but  the  effect  may  be  important  for 
other  species  (eg.  corals)  and  in  a  truly  general  treatment  must  be  taken  into 
account. 

2.  The  equations  derived  here  assume  that  the  flow  forces  experienced  by  an 
organism  are  not  influenced  by  other  organisms  nearby,  by  surface  rugosity,  or 
by  the  presence  of  a  boundary  layer.  These  assumptions  are  probably  valid  for 
barnacles  in  the  high  intertidal  zone  and  a  limpets.  However,  for  other  species 
in  other  habitat  this  assumption  is  clearly  violated,  mussels  being  a  good 
example.  The  mussel  tenacities  measured  in  this  study  were  taken  from  solitary 
mussels,  which  are  occasionally  found  on  Tatoosh  Island,  and  thus  are 
appropriate  for  the  calculations  made  above.  However,  the  vast  majority  of 
mussels  on  Pacific  Northwest  shores  do  not  occur  as  solitary  individuals,  but 
instead  form  large,  ti^itly-packed  beds  (Paine  and  Levin  1981).  Althou^i 
hydrodynamic  forces  may  be  important  in  limiting  the  sizes  of  organisms  in 
aggregations,  as  will  be  discussed  below,  these  forces  are  not  appropriately 
described  by  the  model  presented  here. 

3*  The  model  presented  here  assumes  that  the  coefficients  C^,  C^,  and  C-^ 
are  constant  for  a  given  body  shape.  In  at  least  the  case  of  C^,  the  drag 
coefficient,  this  is  a  gross  oversimplification;  clearly  varies  as  a  function 
of  Reynolds  number,  Re,  (Figure  1)  and,  for  a  given  set  of  flow  parameters,  is 
thus  itself  a  function  of  size.  The  dependence  of  the  inertia  coefficient,  C^, 
and  the  lift  coefficient,  Cp  on  Re  is  less  clear,  as  is  the  dependence  all  of 


these  coefficients  on  period  parameter,  K. 

While  insufficient  data  exist  to  precisely  define  the  manner  in  which  ary 
size-dependent  variation  in  C^,  C^,  and  Cj  affect  the  model,  a  general 
examination  is  possible.  As  the  Reynolds  number  rises  through  the  range  1C)4  - 
1(£  in  steady  state  flow,  C<j  generally  decreases  gradually  (Hoerner  1965,  Vogel 
1981).  ?or  example,  the  drag  coefficient  of  a  smooth  sphere  decreases  from 
approximately  0.47  to  0.10  as  Lj.  increases  from  0.005  m  to  0.025  m  in  a  steady 
flow  velocity  of  10  m/s  (Hoerner  1965).  However,  if  the  surface  of  the  object  is 
rough  or  the  shape  less  symmetrical,  this  effect  is  less  pronounced.  ?or 
example,  a  cylinder  with  surface  roughness  elements  of  height  only  1#  of  the 
diameter  has  a  decreasing  only  from  1 .2  to  approximately  0-9  -  1.0  in  the  Re 
range  10^  -  10^.  further,  as  mentioned  above,  the  drag  coefficient  of  a  plate¬ 
like  shape  decreases  only  slightly  through  a  very  large  range  of  Re.  The  most 
smooth  and  regularly  shaped  organisms  dealt  with  here,  limpets,  have  a 
decreasing  only  from  approximately  0.7  to  0.5  over  a  tenfold  increase  in  Re  (10^ 
to  10  5).  Thus,  while  is  undoubtedly  a  function  of  size,  its  variation  with 
size  for  these  organisms  is  likely  to  be  minimal.  The  variation  in  with  K 
has  not  been  studied  for  biological  objects.  However,  Sarpkaya  and  Tuter  (1974) 
have  shown  that  at  Re  of  5.1 0*  to  105,  the  of  a  cylinder  is  lower  in 

harmonically  oscillating  flow  than  in  steady  flow.  Thus,  the  values  used  here 
probably  overestimate  the  actual  values  at  high  Rs. 

The  only  relevant  C1  va.  Re  data  available  for  shapes  near  surfaces  are 
those  presented  in  Figure  1.  For  the  Re  range  measured,  C]_  appears  to  be 
independent  of  Re,  and  therefore  of  size.  While  these  data  are  far  too 
preliminary  to  draw  any  general  conclusions  they  do  suggest  that  gross 
fluctuations  in  in  the  Re  range  of  interest  are,  unlikely.  The  variation  in 
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with  K  for  biological  objects  is  unknown. 

The  possibility  of  variation  in  as  a  function  of  size  is  difficult  to 
examine.  In  "ideal",  (i.e.  inviscid)  flow,  theory  proclaims  that  is  a  constant 
for  a  given  shape,  regardless  of  size  (Batchelor  1967).  The  method  by  which 
values  were  measured  here  minimizes  the  "non-ideal",  aspects  of  the  flow  by 
measuring  at  very  low  velocities;  and  the  measurements  are  thus  close  to 
those  predicted  by  theory  (Daniel,  1982).  The  nature  of  the  acceleration 
reaction  in  real,  viscid  fluids  is  one  of  the  current  active  areas  of 
investigation  in  fluid  dynamics.  However,  Sarpkaya  and  Tuter  (1974)  have  shown 
for  spheres  and  cylinders  that  for  K  greater  than  approximately  20  (the  range  of 
interest  here)  Cjj,  does  not  vary  substantially.  Sarpkaya  and  Tuter  (1974)  also 
show  that  ^  increases  as  the  Re  increases  in  the  range  from  Re  =  10^  -  10^.  If 
these  results  hold  true  for  biological  objects  near  a  solid  boundary  it  is  quite 
possible  that  values  higher  than  those  used  in  the  above  calculations  are 
present.  Further  research  will  be  necessary  to  resolve  this  possibility. 

We  speculate,  then,  that  when  a  full  knowledge  of  the  variation  in  C<j,C^, 
and  is  available  the  values  used  here  will  be  overestimates,  the  values 
approximately  correct,  and  the  values,  underestimates.  The  net  result  of  a 
model  corrected  for  the  size-dependent  variation  in  coefficients  would  be  in 
each  case  to:  1)  lower  the  estimate  of  the  size  independent  probability  of 
di s lodgement  and  2)  increase  the  L  dP^/di.  Thus,  if  our  speculation  about  the 
size-dependence  of  C^,  and  proves  true,  our  argument  for  the  role  of 
water  acceleration  in  size  limitation  will  be  strengthened  by  the  inclusion  of 
variable  coefficients. 

4.  The  calculation  of  L  dP^/di  depends  in  part  on  the  distribution  of 
breaking  stresses  present  in  a  population.  This  distribution  was  estimated  here 


using  actual  breaking  stress  measurements  which  themselves  contain  some 
experimental  error.  Unless  these  experimental  errors  are  strongly  biased  (e.g. 
the  basal  area  of  small  barnacles  is  overestimated  while  that  of  large  barnacles 
is  underestimated),  and  all  reasonable  precautions  were  taken  to  preclude  such 
biases,  the  experimental  errors  serve  to  increase  the  apparent  variance  of  the 
breaking  stress  distribution.  Prom  an  examination  of  equation  40  and  41  it  can 

be  seen  that  an  increase  in  variance  decreases  the  calculated  L  dP^/dl-  Thus, 
any  future  improvement  in  measuring  technique,  by  giving  a  better  estimate  of 

the  true  breaking  stress  distribution  will,  by  decreasing  the  variance,  increase 
the  estimate  of  L  dP^/dL  for  a  species  and  thereby  strengthen  the  argument  made 
here  for  the  role  of  acceleration  in  size  limitation. 

5-  The  model  proposed  here  assumes  that  adhesive  tenacity  is  independent  of 
the  rate  at  which  a  dislodging  force  is  applied.  This  need  not  be  the  case,  and 
in  general  the  strength  of  biological  materials  increases  as  the  rate  at  which 
they  are  deformed  increases  (Wainwrigit  et^  al^  1976).  However,  the  only  data 
available  concerning  this  point  for  wave-swept  or^nisms  are  those  of  Grenon  and 
Walker  (1982)  who  showed  that  the  adhesive  tenacity  of  the  limpet  Patella 
vulgata  decreases  as  the  dislodging  deformation  is  applied  faster.  Further 
research  is  necessary  before  the  full  effect  on  this  analysis  of  the  rate  of 
deformation  will  be  known. 


Mechanisms  that  Permit  Wave-Swept  Organisms  to  Attain  Large  Size 

Our  predictions  of  mechanical  limits  to  size  of  wave-swept  organisms  have 
been  made  for  solitary,  rigid  animals.  There  are  a  number  of  plants  and  animals 
on  surf-beaten  shores  with  features  that  permit  them  to  attain  sizes  larger  than 
our  model  would  predict.  Among  the  mechanisms  permitting  large  size  are 
aggregation,  posturing,  flexibility,  pruning,  and  reproduction  before  seasonal 


or  rare  large  waves  hit. 


Before  exploring  these  mechanisms  permitting  large  size,  we  should  ask 
whether  the  largest  wave-swept  organisms  (seaweeds  and  corals)  simply  use  a 
stronger  basal  adhesive  and  are  made  of  stronger  material  (i.e.  have  a  higher 
°^naY)  than  their  smaller  neighbors.  Values  for  °uiax  ^or  var>ious  seaweeds 
(e.g.,  Delf  1932,  Koehl  and  Wainwright  1977,  Koehl  1979),  and  corals  (e.g., 
Wainwri^it  et  aL.  1976,  Vosburgh  1977,  1982,  Chamberlain  1978,  this  study)  fall 
within  or  below  the  a^SlX  values  reported  here  for  smaller  organisms. 

Furthermore,  once  the  adhesive  strength  of  the  organism  exceeds  that  of  the  rock 
or  other  organisms  to  which  it  is  stuck,  the  organism  cannot  effectively 
increase  its  ability  to  stay  attached  to  the  shore  by  a  further  increase  in 
strength.  When  dislodging  plants  or  animals,  it  was  found  that  in  many  cases 
the  rock  to  which  the  organisms  were  attached  failed  before  the  organism  came 
unstuck.  For  example,  of  351  R_  glandula  dislodged  at  Mukkaw  Bay,  47$  came  off 
because  the  rock,  rather  than  the  barnacle,  failed. 

As  mentioned  above,  many  intertidal  organisms  occur  in  dense  aggregations. 
Monospecific  stands  may  be  due  to  asexual  reproduction  (e.g.,  clonal  sea 
anemones,  Anthopleura  elegant issima)  or  to  patterns  of  settlement,  growth,  and 
survival  (e.g.,  stalked  barnacles,  Pollicipes  polymerus,  mussels,  M. 
calif omianus,  seaweeds,  palmaefornis).  Furthermore,  on  many  wave-swept 
shores  space  is  a  limited  resource  (see  for  example,  Connell  1961b,  1978,  Paine 
1966,  Dayton  1971,  Quinn  1979,  Paine  and  Levin,  1981)  and  on  such  sites  the 
organism  surrounded  by  bare  substratum  rather  than  a  host  of  other  organisms  is 
rare  indeed.  Because  organisms  in  aggregations  can  be  shielded  by  their 
neighbors,  they  can  be  in  microhabitats  exposed  to  much  lower  flow  velocities 
and  accelerations  than  those  characterizing  mainstream  flow  at  a  site 
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(Wainwright  and  Koehl  1976;  Koehl  1977a,  1977b).  The  water  flow  within  mussel 
beds,  even  on  exposed  coasts,  can  be  so  slow  that  sediment  accumulates  and 
delicate  organisms  can  live  there  (Connell  1972,  Suchanek  1979). 

At  least  in  some  cases  hydrodynamic  forces  are  important  in  removing 
aggregated  organisms  from  the  shore,  large  organisms  or  clumps  being  more 
susceptible  than  small  ones  to  being  blown  away  (e.g.,  Harger  and  Landenberger 
1 971 ,  Connell  1 972,  Paine  and  Levin  1 981 ).  For  our  model  to  be  useful  in 
predicting  the  probability  of  washing  away  of  gregarious  animals  or  plants,  flow 
and  tenacity  measurements  for  organisms  within  a  clump  are  required.  For  some 
organisms  (eg.  Pollicipes  polymerus),  a  more  useful  approach  might  be  to  apply 
the  model  to  the  entire  clump  rather  than  to  individuals  within  it. 

While  considering  possible  mechanical  limits  to  size,  note  that  if  an 
or^mism  on  a  crowded  shore  is  bigger  than  its  nei^ibors,  it  can  stick  out  into 
much  faster  flow  (e.g.,  Wainwright  and  Koehl  1976).  The  difference  in  magnitude 
between  the  mechanical  stresses  experienced  by  the  large  and  small  organisms  in 
this  case  would  be  even  greater  than  our  predictions,  which  were  made  assuming 
that  the  flows  encountered  by  large  and  small  creatures  are  the  same. 

Some  animals  can  actively  change  their  body  shape  or  orientation  in 
response  to  water  flow.  For  example,  the  surge-channel  sea  anemones  Anthopleura 
xanthogrammica  can  actively  adjust  their  height  so  that  they  do  not  protrude 
above  their  neighbors,  and  can  assume  a  pancake-like  posture  (thereby  minimizing 
0^  and  C^)  when  in  rapid  flow  (Koehl  1977a). 

Asymmetric  animals  such  as  mussels  may  grow  with  their  long  axis  oriented 
in  the  direction  of  flow,  thus  reducing  the  added  mass  coefficient.  Orientation 
is  only  useful  when  flow  direction  is  predominantly  along  one  axis,  such  as  in 
surge  channels.  Observations  on  Tatoosh  suggest  that  mussels  are,  in  fact, 
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oriented  with  their  long  axis  parallel  to  flow  in  surge  channels. 

Mobile  animals,  such  as  snails,  can  move  into  protected  microhabitats  on 
wave-beaten  shores.  As  mentioned  above, the  size  of  such  refuges  may  in  some 
cases  impose  an  upper  limit  on  the  sizes  these  animals  can  attain  (e.g.,  Kohn 
1971,  Emson  and  Faller-Fritsch  1976). 

Flexible  organisms  such  as  seaweeds  on  wave-swept  shores  can  attain 
relatively  large  size.  For  example,  macroalgae  such  as  Lessoniopsis  littoralis  in  the 
low  intertidal  zone  on  Tatoosh  Island  can  reach  lengths  (stipe  plus  blades) 
greater  than  1  m  (see  Abbott  and  Hollenberger  1976).  Flexibility  can  decrease 
the  hydrodynamic  forces  organisms  must  resist,  and  hence  can  make  larger  bodies 
mechanically  possible,  by  several  mechanisms:  1)  Tall,  flexible  organisms  tend 
to  be  bent  over  by  moving  water  such  that  their  long  axes  are  parallel  to  the 
flow.  Furthermore,  the  branches  or  lobes  of  flexible  organisms  may  be  stacked 
or  folded  by  moving  water  into  more  compact  shapes  with  lower  drag  coefficients 
(Cfl)  (Fraser  1 962;  Charters  et.al.  1 969;  Koehl  1 977a;  Vogel,  in  press).  Since 
the  drag  force  on  an  organism  depends  on  its  projected  area  normal  to  the  flow 
(Ag)  and  on  its  (eq.  17),  the  drag  on  a  flexible  organism  can  be  lower  than 
on  a  rigid  organism  of  the  same  size  and  original  shape.  The  bending  over  or 
deforming  of  flexible  organisms  in  waves  can  also  reduce  the  acceleration 
reaction  force  (eq.  19);  the  C^s  of  flattened  bodies  parallel  to  the 
acceleration  is  lower  than  those  of  more  spherical  bodies,  which  in  turn  are 
lower  than  the  2m  *9  of  flattened  bodies  normal  to  the  acceleration  (Daniel 
1983).  Furthermore,  if  a  tall  flexible  organism  is  surrounded  by  other 
organisms  or  is  thin  (i.e.  very  flat  or  narrow)  compared  with  the  boundary  layer 
(  S  ,  eq.  16)  along  the  substratum,  it  can  be  bent  over  by  moving  water  into  a 
position  closer  to  the  substratum  where  it  encounters  slower  flow.  2)  Flexible 


organisms  that  are  long  enough  can  avoid  bearing  hydrodynamic  forces  at  times 
during  a  wave  when  they  are  likely  to  be  greatest  (i.e.  when  du/db  and  u  are 
highest,  eqs.  19  and  17).  For  example,  a  very  floppy,  long  alga  such  as 
Durvillea  antarctica  moves  with  the  water  in  a  wave  as  it  rushes  up  the  shore; 
the  thallus  of  the  plant  is  not  substantially  stretched  until  the  alga  is  strung 
out  in  the  direction  of  flow  and  is  no  longer  free  to  move  with  the  water.  If 
such  a  flexible  organism  is  long  enough,  the  water  may  be  slowing  down  or  have 
reversed  direction  before  the  slack  in  the  plant  can  be  taken  up.  Thus,  for 
long,  floppy  or^uiisms  in  oscillating  flow,  an  increase  in  length  can  lead  to  a 
decrease  in  the  forces  the  organism  has  to  bear.  3)  Deformable  organisms  are 
good  "shock  absorbers".  A  hydrodynamic  force  on  such  an  organism  bends  or 
stretches  the  body,  and  not  until  the  body  is  fully  deformed  is  the  entire  load 
transmitted  to  the  area  of  attachment  to  the  substratum.  Thus,  hi$i  forces  of 
short  duration  (such  as  occur  an  organisms  in  waves)  can  be  damped  out  by 
bendable,  soft,  or  stretchable  bodies  (see  Pain  1968;  Thomson  1981). 

That  flexibility  can  be  an  effective  mechanism  for  withstanding  wave  action 
was  clearly  illustrated  when  Hurricane  Allen  hit  the  Jamaican  coral  reefs: 
rigid,  brittle  organisms  such  as  stony  cora1s  suffered  more  breakage  than  did 
flexible,  deformable  organisms  like  gorgonians  (Woodley  et  al.  1981). 
Nonetheless,  we  should  also  point  put  that  in  certain  cases,  flexiblity  may 
reduce  the  likelihood  that  an  organism  will  survive  in  a  wave-beaten  habitat. 

If  a  flexible  organism  flaps  (like  a  flag)  in  the  flow,  it  can  experience  higier 
drag  than  a  rigid  structure  of  the  same  size  and  shape  (e.g.  Hoerner  1965). 
Furthermore,  if  structures  are  flexible  enough  to  oscillate  in  waves,  they  can 
perish  by  fatigue  fracture  (e.g.  Weigel  1972).  The  behavior  of  flexible 
organisms  in  waves  is  a  complex  mechanical  problem  that  merits  further  analysis. 


The  heaviest  wave  action  generally  accompanies  seasonal  storms  (e.g. 
Schwenke  1971,  Harger  1972,  Menge  1972,  Paine  and  Levin  1981)  or  isolated  events 
like  hurricanes  and  typhoons  (e.g.  Stoddart  1962,  1963,  Endean  1976,  Randall  and 
Eldridge  1977,  Smith  and  Harrison  1977,  Vosburgh  1977,  Woodley,  et  al.,  1981, 
Hi^ismith  1982).  Some  sessile  organisms  (e.g.,  certain  algae  and  colonial 
invertebrates)  that  grow  "too  large",  during  non-stormy  intervals  are  pruned  to  a 
smaller  size  during  periods  of  heavy  wave  action  rather  than  being  completely 
ripped  off  the  substratum.  The  places  at  which  pieces  of  the  organism  break  off 
may  be  determined  by  the  construction  of  the  organism  itself  (e.g.  Highsmith 
1982)  or  by  the  activities  of  grazers  on  and  bioeroders  in  the  organism  (e.g., 
Black  1 976:  Santelices  et  al.  1 980,  Highsmith  1 982).  The  "stump",  that  is  left 
behind  after  such  pruning  can  continue  to  live  and  may  grow  to  large  size  again 
during  subsequent  intervals  of  calm  water.  In  some  cases,  not  only  can  the 
stump  survive,  but  so  can  broken-off  fragments.  "Programmed  fragmentation",  of 
this  sort  has  been  proposed  ’as  a  mechanism  of  assexual  reproduction  and 
dispersal  for  many  corals  (e.g.,  Chamberlain  1978,  Tunnicliffe  1980,  1981, 
Highsmith  1982). 

In  habitats  where  large  waves  occur  only  seasonally  or  during  rare  storms, 
organisms  that  grow  rapidly  during  calm  intervals  can  become  too  large  to 
survive  the  next  onslaught  of  violent  water  movement.  However,  if  such 
organisms  can  produce  gametes  or  spores  before  they  are  blown  away,  and  if  their 
propagules  can  recolonize  the  habitat,  then  selection  for  smaller  size  by  wave- 
induced  mortality  should  not  be  strong.  An  example  of  such  a  rapidly-growing, 
early-reproducing  large  organism  that  thrives  on  exposed  shores  that  have 
predictable  seasonal  periods  of  heavy  wave  action  is  the  seaweed  Durvillea 
antarctica  (Santilices  et  al.  1980). 
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To  this  point  we  have  dealt  solely  with  mechanical  limits  to  size,  and  have 
treated  large  size  as  if  it  were  an  unquestioned  "good  thing?',  i.e.  advantageous 
to  the  organism.  This  is  not  necessarily  so,  and  in  order  to  place  the  ideas 
proposed  here  in  a  context  that  will  be  useful  in  the  examination  of  real  wave- 
swept  communities,  we  briefly  examine  the  biological  advantages  and 
disadvantages  of  size. 

Is  Big  Good? 

Large  size  may  be  advantageous  to  an  organism  by  increasing  its  persistence 
on  the  shore,  by  increasing  its  reproductive  output,  or  both.  A  number  of 
advantages  to  large  size  have  been  suggested: 

large  individuals  are  better  able  to  ward  off  some  types  of  pedators  than 

are  small  individuals  (e.g.  Gould  1966,  Connell  1970a,  Stanley  1973,  Zaret  and 

Kerfoot  1975).  A  number  of  cases  of  benthic  or^nisms  escaping  in  size  from 

» 

predation  have  been  documented  (Ansell  I960,  Paine  1965,  1976a,  Dayton  1971, 
Palmer  1979,  and  a  number  of  other  studies  reviewed  by  Connell  1972  and  by 
Vermeij  1978).  Furthermore,  in  some  cases  inferior  competitors  can  escape  in 
size  from  being  overgrown  by  members  of  other  species  (Sebens  1980).  Tall 
solitary  organisms  are  often  less  susceptible  than  smaller  ones  to  being 
overgrown  by  encrusting  colonial  animals  (Jackson  1977). 

Because  large  organisms  have  a  lower  surface  area  to  volume  ratio  than  do 
smaller  creatures,  the  large  organisms  .are  better  insulated  from  fluctuations  in 
external  conditions  and  can  regulate  their  internal  environment  more  easily 
(Gould  1966,  Stanley  1973).  Many  organisms  exposed  to  wave  action  are  intertidal 
and  thus  are  often  exposed  to  drying  in  air  and/or  heating  in  the  sun.  That 
large  intertidal  organisms  survive  exposure  to  air  longer  than  do  smaller 


members  of  the  same  species  has  been  demonstrated  (e.g.  Johnson  and  Shick  1977). 
Although  sessile  individuals  must  survive  being  small  in  the  same  spot  that  they 
occupy  when  large  (Paine  1981),  daytime  low  spring  tides  during  warm  sunqy 
weather  are  seasonal  occurrences  at  many  sites  -  an  organism  at  such  a  site 
could  conceivably  settle  and  grow  beyond  desiccation-vulnerable  small  size 
before  the  harsh  tides  occur. 

large  organisms  may  have  a  greater  likelihood  of  surviving  damage  than 
small  ones.  Plants  and  animals  on  wave-beaten  shores  are  sometimes  bashed  by 
water-borne  projectiles  such  as  logs  (Payton  1971)  or  coral  rubble  (Endean  1976, 
Randall  and  Eldredge  1977,  Highsmith  et.al  1980).  Although  large  corals  are  more 
likely  to  suffer  some  damage  by  projectiles  than  are  small  ones  (Woodley,  et. 
al.  1 981 ),  large  corals  are  less  likely  than  small  ones  to  be  detached  or 
destroyed  when  hit  by  a  projectile  (Highsmith  1981,  Woodley  et.al.  1981). 
Furthermore,  large  pieces  of  coral  that  are  broken  off  have  a  greater  chance  of 
surviving  than  do  small  ones  (Highsmith  1980,  1982,  Highsmith  et.al.  1980). 
Converse1  v,  of  course,  if  an  animal  is  small  enough  to  live  in  a  crevice,  its 
probability  of  being  scraped  off  the  substratum  by  a  projectile  is  very  low 
(Connell  1972,  Paine  1981). 

As  mentioned  previously,  size  may  directly  affect  repoductive  output.  Largs 
animals  have  more  "room",  (volume)  in  their  bodies  in  which  gonadal  tissue  might 
develop  (Sebens  1982).  Gould  (1966)  has  suggested  that  large  invertebrates  can 
have  more  progeny  per  brood  than  small  ones,  and  Menge  (1974)  has  demonstrated 
that  the  fecundity  of  large  Leptasterias  (starfish)  is  greater  than  that  of 
small  ones.  This  relationship  is  quite  general  among  invertebrates;  for  a  review 
see  Giese  and  Pearse  (1974). 

Another  advantage  of  increasing  body  size  is  that  new  catagories  of  food 
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can  become  available  to  an  organism  if  it  becomes  big  enough  (e.g.  Gould  1966, 
Stanley  1973,  Wilson,  1975).  Examples  of  wave-swept  organisms  that  can,  once 
they  attain  a  large  size,  capture  and  ingest  new  types  of  food  are  provided  ty 
Menge  (1972)  and  Sebens  (1979,  1981,  1982).  As  Armstrong  (1977)  has  pointed  out, 
large  organisms  that  can  take  big  food  items  may  have  a  competitive  advantage 
over  small  ones  when  the  food  resource  distribution  is  skewed  towards  large 
size. 

large  or^nisms  are  more  efficient  metabolically  (i.e.  have  a  lower 
metabolic  rate  per  unit  body  mass)  than  small  organisms  (e.g.  Gould  1966, 
Schmidt-Nielsen  1974).  Furthermore,  big  mussels  can  achieve  their  maximum  growth 
rate  on  a  lower  relative  food  ration  than  can  small  ones  (Griffiths  and  King 
1979). 


Biological  Factors  Limiting  the  Size  of  Wave-Swept  Organisms 

In  spite  of  the  many  benefits  of  being  big,  a  number  of  biological  factors 
can  counteract  the  selective  advantage  of  large  size  or  can  locally  prevent 
individuals  from  attaining  the  maximum  size  that  is  both  biologically  and 
mechanically  possible  for  their  species.  For  example,  some  snail  and  starfish 
predators  have  been  found  to  preferentially  eat  large  rather  than  small 
individuals  of  some  barnacle  species  (Connell  1961a,  1970b,  Paine  1966,  1981, 
Palmer  1980),  and  may  thus  remove  the  big  individuals  from  a  population.  If  a 
particular  type  of  organism  must  hide  (from  dessication,  waves,  or  predators, 
for  example),  the  sizes  of  refuges  available  at  a  particular  site  may  limit  the 
maximum  sizes  of  the  organisms  (e.g.,  Kohn  1971,  Emson  and  Faller-Fritsch  1976, 
Paine  1981);  organisms  too  big  to  fit  into  the  refuges  perish.  The  allocation 
of  resources  to  the  repair  of  wave-inflicted  damage  ratner  than  to  growth  may 


also  limit  the  size  of  wave-beaten  creatures  (e.g.,  Ebert  1968). 

The  size,  and  thereby  the  reproductive  output,  an  organism  attains  also 
depends  in  part  upon  the  quantity  and  quality  of  food  available  to  it  in  a 
particular  habitat  relative  to  its  metabolic  demands  (e.g.,  sea  anemones:  Sebens 
1982;  gastropods:  Prank,  1965,  Paine  1965,  1969,  Sutherland  1970;  bivalves: 
Stanley,  1975,  Griffiths  and  King  1979;  starfish:  Smith  1940,  Paine,  1976a;  sea 
urchins:  Ebert  1968,  Vadas  1977;  fish:  Kerr  1971).  Many  intertidal  animals  may 
have  to  stop  feeding  when  exposed  to  air  or  when  wave  action  is  severe,  and  thus 
may  spend  a  smaller  fraction  of  each  day  eating  than  subtidal  animals.  Reduced 
feeding  time,  in  combination,  in  some  cases,  with  lower  food  availability  or 
increased  metabolic  rate  (due  to  warming  up  when  out  of  the  water),  have  been 
cited  as  responsible  for  the  observed  trend  that  high  intertidal  animals  are 
smaller  than  low  intertidal  ones  (Seed  1969,  Sutherland  1970,  Connell  1961a, 
1972,  Dayton  1971,  Vermeij,  1972,  1978,  Paine  1974,  Suchanek,  1979,  Sebens 
1980).  This  pattern  of  size  distribution  has  also  been  ascribed  to  the 
migration  of  larger  animals  to  lower  positions  on  the  shore  (Paine  1969, 

Bertness  1977).  Note  that  a  few  intertidal  species  show  the  opposite  trend  in 
size  distribution  (i.e.,  larger  individuals  higher  on  the  shore)(e.g.,  Prank 
1965).  Various  special  mechanisms  responsible  for  thet  3  exceptions  are  listed  by 
Vermeij  (1978). 

Although  large  animals  have  a  lower  metabolic  rate  per  unit  mass  than  small 
canes,  their  total  metabolic  rate  per  individual  is  greater.  If  the  metabolic 
rate  of  an  animal  increases  as  a  greater  power  of  body  mass  than  does  the  rate 
at  which  an  animal  can  take  in  energy,  the  scope  for  growth  ’(sensu  Sebens,  1979) 
rises  to  some  maximum  and  then  decreases  as  body  size  increases  (e.g.,  mussels: 
Vahl  1973,  Thompson  and  Bayne  1974,  Griffiths  and  King  1979;  cnidarians:  Sebens 
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1979,  1981,  1982).  For  such  organisms  an  upper  size  limit  exists  at  the  point 
where  scope  for  growth  equals  zero.  Sebens  (1979,  1981)  suggests  that  such 
oganisms  should  get  no  larger  than  the  "optimal  size",  at  which  scope  for  growth, 
and  hence  gonad  production,  is  maximized.  This  optimal  size  is  smaller  than  the 
maximal  size,  and  is  a  function  of  habitat  quality* 

Applications  of  the  Model  of  Size-Dependent  Mechanical  Failure 

In  light  of  these  many  concurrently  operating  factors  affecting  body  sizes 
in  a  population  how  can  the  proposed  here  be  usefully  applied? 

Our  analysis  provides  a  quantitative  method  for  deciding  whether  wave 
forces  alone  may  limit  the  size  of  a  particular  species  at  a  particular  site. 
From  the  estimated  or  measured  values  for  the  parameters  of  the  model,  the 
probability  of  dislodgement  and  the  size-specific  increment  in  probability  of 
dislodgement  can  be  calculated.  If  these  values  are  high,  it  is  possible  that 
the  size  of  that  species  is  limited  at  that  site  by  the  wave  forces  encountered, 
and  further,  experimental  efforts  may  be  taken  to  confirm  or  refute  thi3 
conclusion.  In  such  a  case  where  P^  and  I*  dP^/dl  are  high,  an  examination 
solely  of  the  possible  biological  factors  limiting  size  would  be  ill-advised. 
Conversely  if  and  L  dP^/dl  are  low,  wave  forces  may  reasonably  be  discounted 
and  further  efforts  focused  on  the  role  of  biological  interactions. 

Even  in  those  cases  where  wave  forces  may  be  deemed  to  be  only  marginally 
involved  in  determining  the  upper  size  limit,  a  consideration  of  these  forces 
mgy  prove  useful  in  examining  the  optimal  (as  opposed  to  the  maximal)  size  that 
a  species  attains.  The  hypothesis  proposed  here  provides  a  mechanism  whereby 
larger  organisms  are  exposed  to  a  greater  probability  of  death.  This  size- 
specific  mortality  could  conceivably  serve  to  decrease  the  size  at  which  an 


organism  optimizes  its  intrinsic  rate  of  increase  (sensu  Roff  1931)  and  could 
thereby  affect  selection  for  smaller  individuals. 

Alexander  (1981)  has  outlined  a  theory  whereby  the  safety  factors  with 
which  biological  structures  are  constructed  can  be  predicted.  Data  collected 
from  examination  of  animal  structure  in  regards  to  optimal  size  may  also  be  used 
to  predict  optimal  safety  factors. 

The  methods  outlined  here  may  prove  useful  in  examining  the  morphological 
diversity  of  organisms  present  in  a  given  habitat.  Vermeij  (1978)  has  suggested 
that  the  greater  diversity  of  snail  shell  types  found  in  protected  as  compared 
with  exposed  habitats  is  due  to  the  limits  on  shell  form  imposed  by  wave  action. 
Similarly,  we  mi^it  expect  to  find  a  greater  diversity  of  form  among  small  than 
among  large  organisms  on  wave-swept  shores.  Stanley  (1973),  Bonner  and  Horn 
(1982),  and  Horn  et  al.  (1982)  have  mentioned  various  other  physical  constraints 
on  form  that  affect  the  morphological  diversity  of  large,  but  not  small, 
or^misms.  These  propositions  can  now  be  quantitatively  tested. 

The  ideas  proposed  here,  by  quantifying  the  "exposure",  (ie.  the  probability 
of  dislodgement)  of  organisms,  may  prove  useful  in  several  aspects  of  behavioral 
and  community  ecology.  Because  wave  action  can  limit  the  time  a  motile  organism 
can  forage  without  being  in  danger  of  washing  away  (e.g.,  Menge  1972),  and 
because  large  animals  (both  predators  and  prey)  are  more  likely  to  be  washed 
away  than  small  ones,  the  foraging  strategies  of  large  vs.  small  animals  on 
protected  vs.  exposed  shores  might  be  expected  to  be  different.  The  role  of 
disturbance  in  maintaining  the  diversity  of  communities  has  been  discussed  by  a 
number  of  authors  (e.g.,  Levin  and  Paine  1974,  1975,  Paine  and  Levin  1981, 
Connell  and  Slayter  1977,  Connell  1978,  Sousa  1979a,  1979b,  Quinn  1979,  and  many 
others  cited  therein).  Waves,  which  are  an  important  agent  of  disturbance  on 
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many  rocky  shores  and  coral  reefs,  are  more  likely  to  remove  large  than  small 
organisms  from  the  substratum.  Such  differential  suscepti  ility  of  organisms  to 
disturbance  may  have  important  consequences  to  the  age  structure  and  species 
composition  of  wave-swept  communities.  Purthermore,  the  time  required  for  the 
particular  primary  space-holders  in  a  community  to  grow  to  wave- vulnerable  size, 
coupled  with  the  temporal  pattern  of  periods  of  heavy  wave  action  at  a  site, 
should  be  of  basic  importance  in  the  dynamics  of  the  community  (e.g.,  Paine  and 
Levin  1 981 ). 

Size-dependent  predation  can  have  important  effects  on  population  and 
community  structure  (e.g.,  Brooks  and  Dodson  1965,  Galbraith  1967,  Zaret  and 
Kerfoot  1975).  Wave  action  that  selectively  removes  large  individuals  from  a 
population  should  have  similar  effects.  Note  that  the  more  exposed  a  habitat, 
the  smaller  the  size  at  which  it  is  highly  probable  that  a  particular  type  of 
organism  will  be  washed  away. 

Sessile  organisms  often  provide  a  habitat  in  which  other  organisms  live. 
Often  the  larger  the  sessile  organisms,  the  greater  the  diversity  of  the 
community  of  organisms  living  amongst  them  (e.g.,  MacArthur  and  MacArthur  1961, 
Smith  1972,  Paine  1976a,  Suchanek  1979)*  Thus,  the  sizes  at  which  sessile 
organisms  are  likely  to  wash  away  in  habitats  of  different  wave  exposure  should 
have  important  effects  on  overall  community  structure. 


SUMMARY 


1.  There  are  many  advantages  to  large  size. 

2.  Despite  these  advantages  most  wave-swept  organisms  are  small. 

5-  The  accele rational  component  of  flow  in  breaking  waves  places  a  size- 
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dependent  mechanical  stress  on  wave-swept  organisms.  This  size-dependent  stress 
is  hypothesized  to  be  a  limiting  factor  in  the  size  to  which  many  wave-swept 
organisms  grow. 

4.  This  hypothesis  is  tested  using  various  intertidal  species,  and  is 
conjectured  to  be  valid  for  the  limpets,  Ct_  pelta  and  1L_  scutum,  the  urchin,  S. 
purpuratus,  the  mussel,  californianus  (when  solitary),  and  the  hydrocoral,  M. 
complanata.  For  the  snails  discussed  (T^  emarginata,  T^_  canaliculata,  and  L. 
scutulata)  the  mechanical  environment  poses  an  indirect  constrain  on  body  size 
and  feeding  behavior.  For  the  acorn  barnacle,  S^  cariosus,  mechanical  limits 
appear  to  be  less  of  a  factor  in  limiting  size. 

5.  A  definitve  answer  regarding  the  possibility  of  mechanical  limitations 
to  size  depends  on  factors  that  remain  to  be  measured:  long-term  wave  force 
probabilities,  and  life  history  parameters. 

6.  Flexibility,  posturing,  pruning,  and  life  cycle  strategies  provide 
possible  mechanisms  for  circumventing  the  sizo  limitations  due  to  acceleration. 

7.  The  hypothesis  proposed  here  may  prove  useful  in  examining  the  ecology 
of  wave-swept  communities  by  providing  a  quantitative  means  for  measuring  the 
importance  of  mechanical  factors  in  the  functioning  of  organisms. 
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APPENDIX  1 


Symbol 

Definition 

Units 

Eq.  where  first 

*b 

Basal  area  of  model  or  ganism 

m^ 

20 

he 

Critical  area 

m2 

2 

4 

Projected  area  (direction  of  flow) 

m2 

1 

iss 

Maximum  projected  area 

m2 

42 

Ap.min 

Minimum  projected  area 

m2 

42 

B 

Proportionality  constant 

— 

3 

5d 

Drag  coefficient 

— 

7 

9l 

Lift  coefficient 

— 

8 

5m 

Added  mass  coefficient 

— 

6 

4 

Center  of  pressure 

— 

11 

D 

Specimen  length  (coral) 

m 

13 

d 

Distance  from  neutral  axis 

m 

23 

? 

Loading  force 

N 

13 

4 

Acceleration  reaction 

N 

4 

7d 

Drag  force 

N 

1 

Zdf 

Total  force  in  direction  of  flow 

N 

4 

4 

Lift  Force 

N 

8 

& 

Acceleration  due  to  gravity 

m/s2 

14 

H 

Wave  height 

m 

14 

Srma 

Root-mean  square  wave  height 

m 

46 

4 

Pressure  head  at  port  i 

m 

11 

i 

Second  moment  of  area 

m^ 

13 

90 


Proportionality  constant 
Period  parameter 
Characteristic  length 
Height  of  coral  blade 
Radius  of  animal  or  model 
Maximum  radius 
Mass 

Moment  (force  x  distance) 

Number  of  waves  encountered 

P,j  due  to  acceleration 

Probability  of  encountering 
acceleration  A 

Size-independent  P^ 

Probability  of  dislodgement 

Cumulative  Pfl 

Probability  that  a  wave  has  a 
height  >  H' 

Probability  of  survival 

Cumulative  probability  of  survival 

Probability  of  encountering 
velocity  U 

Probability  of  encountering 
velocity  in  range 

Probability  that  a  wave  has  a 
velocity  >  U' 

Radial  distance  to  port  i 

Reynolds  number 

Standard  deviation  of  breaking 
stress 


eg, max 
n 

n.max 


Relative  thickness  (of  coral  blade) 

— 

Period  of  oscillation 

s 

Water  velocity 

m/s 

Maximum  water  velocity  in 

oscillatory  flow 

m/s 

Volume  of  displaced  fluid 

m? 

Distance  from  neutral  axis 

m 

Relative  width  (of  coral  blade) 

m 

Distance  from  leading  edge 

m 

Water  depth 

m 

Angle  relative  to  mainstream  velocity 

degrees 

Added  mass  coefficient 

— 

Boundary  layer  height 

m 

Kinematic  viscosity  of  water 

m2/s 

Cone  angLe 

degrees 

Density  of  water 

kg/m? 

Stress  (force/area) 

N/m2 

Average  breaking  stress 

N/m2 

Equivalent  maximum  stress 

N/m2 

Normal  stress 

N/m2 

Maximum  normal  stress 

N/m2 

Maximum  shear  stress 

N/m2 

APPENDIX  2 

Any  change  in  the  constants  or  coefficients  appearing  in  eqs.  30  and  31 
causes  the  estimate  of  maximum  size  to  vary.  The  relative  sensitivity  of  Lj.  ma} 
to  variations  in  C<j,  Cjn,  Cj,  U,  CT  max,  and  dU/dt  may  be  assessed  as  follows: 
Ifr  max  33  defined  by  eq.  26  is 


ijymax  -  (  ^"a.max  ”  ^lpU^)/(^  fr?  dU/dt) 


and  for  Kj  =  [1/3]  C<j  and  £2  =  C8/93  (values  for  a  hemisphere) 


^r,max  -  (^s.max  “  [e/3]  Cd  U2)/[8e/9]  C*  dU/dt 


The  total  change  in  Lj.  max  is 


^r>max  =  f  ^  — r>max^<7s>max  ^  d  c5simax 
+  Id^max/^  2d^  dC^ 

+  ^  ^r,max/  *  1  dCj^ 

+  I?  Iirtmax/  *  £m  J  dC* 

+  ^^r.max/3  dU/dt]  d2U/dt2 

?or  small  errors  such  that 

d^r,max  jjr,max  ^  ^r,max  (a^  an>X  ^s,max  £d>  ^m>  U»  ^d  dU/dt) 
d°s,max  ~A‘C-^s,max  ^  ^s.max  (a^  lir.raax’  £d»  £m»  dU/djt) 


93 


etc.  for  C,3,  C^,  U,  and  dU/dt 

we  may  linearize  the  above  equation  as: 

A  Iirtmax  ^  -^max/^  P's ,  max  3A  Qrs,  max 

+  ^Ilr,max/^  C.d.3  A  £d 
+  [^jrtmax/  3  Cjn3  A  Cjm 
+  kjfr, max/a  U]  a  U 
+  [^ia^max/^dU/dt]  a  dU/dt  . 

This  defines  the  absolute  variation  in  Lr>max;  of  more  interest  is  the  relative 
variation  a  hr,maxAytmax: 

A  ^r,maxA;rtmax  =  k  j^max/d  g~stmax3  g~3>maxA^r,max3 

+^Iir,maxA^d3  [^C^/L^  max  3 

4^jjr,max/^^m3  ,max3 
+bir.max/>gl  jjnax  3 
a  jfrt max/d dU/ d^t]  [ a  dU/ dt, Air ^ max ] 

When  these  partial  derivatives  are  taken  and  A=  5(  P^>max  “t^/3]  33d  U2)  j 

A  33r, max/^r , max  =  ^,max^  +  [U^p  A  C^3/3A  +  [^  £d  U  aU]/3A 
+  [*dU/dt]/[dU/dt]  +  a  C^/Cjn 

Thus  a  1C*  variation  in  dU/djt  or  will  always  produce  a  IQ#  variation  in 
Lp^ay.  In  contrast,  the  relative  variation  in  Lj^max  cause(i  hy  a  10#  variation 
in  crs,max>  £d>  or  31  depends  on  the  magnitude  of  P^may,^,  and  U.  ?or  example 
a  IQ#  variation  in  causes  a  0.14#  variation  in  I^max  ^or  a  wall  streamlined 


or^mism  with  a  very  strong  adhesive  under  severe  flow  conditions  (U  =  20  m/s, 
Cj}  =  0.1,  ^jinax3  =  10^  N/m2).  In  contrast,  a  10$  variation  in  C<j  will  cause  a 
58$  change  in  Lp>max  for  a  bluff  organism  with  a  very  weak  adhesive  under  mild 
flow  conditions  (U  =  5  m/s,  &  =  1.0,  <T  max  ^  N/m2).  A  typical  example  from 
this  study  would  be  the  limpet  digitalis  while  moving  in  moderate  to  hi$i 
flow  (U  =  10  m/s,  Ofl  =  0.525,  <3~s>max  1.29  x  105  N/m2)  where  a  10$  variation  in 
0^  would  result  in  a  1.6$  variation  in  Lr>max» 

A  similar  sensitivity  analysis  can  be  carried  out  for  eq.  27  and  is 
presented  here  without  derivation. 

B=[  32/9*]  e  Cm  dU/dt 

£r,max^r.max  -  *  t0'028  »2e  +  [(4/3ti)eU2^Cj/B] 

+  [2  Ue  (0.028  C!  +  (4/3n)C^)  aU/B] 

+  [  -  l2e  (1/2  -  0.472  +  4/3*)  a  C^/tB  C^2]] 

+  [  (  n  -  U2  (1/2  -  0.472  +  4/3 it))  [AdU/dt]/[B  (dU/dt)2] 
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Adhesive  tenacity  values  and  hydrodynamic  parameters  for  a  variety  of  Wave-swept  organisms 
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TABLE  2.  Shape  constants  for  various  standard  shapes. 


Hemisphere 

Cone  . 

Cylinder 

(1/3)  ^ 

(2/3  )B  ^ 

(4/3  )B  C 

(8/9)  ^ 

(4/9)B  ^  ' 

( V3)B  ^ 

(1/2)  C, 

(1/2)  C, 

d/2)  Cj 

-.472  C, 

-.472 

-.472  C, 

(4/311)  C,, 

(2-(2/[2]1/2)B2Cd 

(2/ir)B2  C 

(32/94)  4, 

(4/3  it)  B2  ^ 

2B2  C 

—  -m 

=  Height/radius 


TABLE  3.  Approximate  maximum  sizes  for  organisms  used  in  this  study 


SPECIES 

Semibalanus  cariosus 
Balanus  glandula 
Balanus  nubilis 
Colli sell a  digitalis 
Colli sell a  pelta 
Notoacmaea  scutum 
Mytilus  californianus 

i 

Strong 1 y 1 ocen t rotu s 
purpuratus 

Millepora  comp! ana ta 
Thais  emarginata 
Thais  canal iculata 
Littorina  scutulata 


DIMENSION 


L 


LENGTH 

^  <  -? 

,  i'  1.85  x  10*m 

■& "  'iP'\ 

1.00  x  I0"2m 
4.75  x  I0“2m 
0.85  x  10-2m 
1.50  x  I0"2m 

I. 50  x  I0”2m 

II. 0  x  lo~2m 


3.70  x  I0~2m 
14.0  x  I0~2m 
2.00  x  I0"2m 
2.00  x  10"2m 
0.50  x  10~2m 
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Figure  1.  Drag  and  life  coefficients  (C^  and  C^)  are  plotted  against  log  (Re) 
for  a  variety  of  organisms  and  standard  shapes:  acorn  barnacle,  cariosus 
(closed  circles);  limpet,  (L_  pelta  (open  circles);  hemisphere  (open 
triangles);  sphere  (closed  triangle).  A  lift  force  occurs  only  if 
hydrostatic  pressure  can  be  transmitted  to  the  basal  surface  of  the 
organism.  Thus  an  acorn  barnacle  glued  to  the  rock  with  a  solid  adhesive 
(as  in  nature)  would  not  experience  a  lift.  The  Cj  values  shown  here  were 
measurd  using  a  barnacle  replica  separated  from  the  wall  by  a  thin  layer  of 
water,  an  unnatural  condition,  and  are  presented  solely  to  show  that  for 
such  shapes  does  not  change  substantially  with  Re. 

Figure  2.  A  schematic  diagram  of  the  experimental  method  for  determining  the 
location  of  the  center  of  pressure.  A  plastic  cone  is  immersed  in  flowing 
water.  The  row  of  ports  is  placed  at  some  angle,  ,  relative  to  the 
direction  of  flow;  rj  is  the  distance  from  the  apical  port  (_Tq)  to  the  first 
port.  A  manometer  is  used  to  measure  the  difference  in  pressure  between  each 
port  and  the  underside  of  the  cone. 

Figure  3-  Adhesive  tenacity  is  plotted  against  shell  length  for  the  mussel  M. 
califomianus.  The  equation  which  describes  these  data  is:  Force  =  1.65  x 
10^  L  1.95  where  Force  is  expressed  in  N  and  L  in  m  (r=0.86). 

Figure  4.  Volume  is  plotted  against  shell  length  for  the  mussel  M. 

califomianus.  The  equation  which  describes  these  data  is:  Volume  =  0.061  L 
2.86  where  Volume  is  expressed  in  m3  and  L  in  m  (r=0.98). 

Figure  5.  Water  velocity  (broken  line)  and  acceleration  (dashed  line)  during  the 
initial  portion  of  a  typical  wave,  and  the  hydrodynamic  force  (solid  line) 
they  impose  on  an  acorn  barnacle  are  plotted  against  time. 
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re  6;_A  model  of  the  hydrodynamic  forces  imposed  on  a  hemisphere  attached  to 
a  planar  surface.  The  total  force  in  the  direction  of  flow,  j^f,  places  a 
shear  force  (F^)  and  a  moment  (F^f.  b)  on  the  hemisphere.  The  lift  force, 
?i ,  imposes  a  moment  (F^f-c)  on  the  hemisphere.  C  is  the  centroid  and  N  is 
the  location  of  the  neutral  axis  (a  line  perpendicular  to  the  plane  of  the 
drawing).  Lp  is  the  radius  of  the  hemisphere.  .  The  components  of  ?df  are 
drag,  Pp|,  and  the  acceleration  reaction,  F^ 

^e  7.  Predictions  for  the  urchin  purpuratus  exposed  to  a  variety  of  water 
velocities  (1,5,10,  and  20  m/s),  (a)  Probability  of  dislodgment  (P^)  is 

plotted  against  water  acceleration  (bp, max  =  ^.037  m).  The  size-specific 
increment  in  risk  (L  dP/dl)  is  plotted  against  size  (radius  of  urchin)  for 
an  animal  exposed  to  (b)  500  m/s^  and  (c)  1000  m/s^.  The  maximum  observed 
size  is  shown  by  the  vertical  dashed  line. 

re  8.  Predictions  for  the  limpets  C.  pelta  (a),  N.  scutum  (b),  and  C. 
digitalis  (c)  exposed  to  a  variety  of  water  velocities  (1  to  10  m/s). 
Predictions  are  based  on  adhesive  tenacities  for  moving  (solid  lines)  and 
statinary  (dashed  lines)  limpets.  Left  panels:  probability  of  dislodgsment 
(Pjj)  is  plotted  against  water  acceleration  for  limpets  of  the  observed 
maximum  size  (Lj.  max).  Righ  panels:  the  size-specific  increment  in  risk  (L 
dPfl/dL)  is  plotted  against  size  for  limpets  exposed  to  a  water  acceleration 
of  500  m/s^.  The  observed  maximum  size  is  shown  by  a  vertical  broken  line. 
;e  9^  Predictions  for  the  mussel,  californianus  exposed  to  variety  of 
water  velocities  (1  to  20  m/s)  oriented  both  broadside  (upper  panels)  and 
end-on  (lower  panels)  relative  to  the  direction  of  flow.  Left  panels:  the 
probability  of  dislodgement  (P^)  is  plotted  against  water  acceleration  for 
mussels  of  dimension  Lmay.  Right  panels:  the  size-specific  increment  in 
risk  is  plotted  against  size  for  mussels  exposed  to  a  water  acceleration  of 
500  m/s^«  Lmax  is  the  maximum  size  of  solitary  mussels  observed  on  Tatoosh 


Island. 


Figure  10.  Predictions  for  the  fire  coral,  complanata  exposed  to  a  variety  of 
water  velocities  (1  to  20  m/s).  Left  panel:  the  probability  of  breaking  is 
plotted  a^iinst  water  acceleration  for  a  coral  of  blade  .length  0.14  m. 

Right  panel:  the  size-specific  increment  in  risk  is  plotted  against  blade 
length  for  a  coral  exposed  to  a  water  acceleration  of  500  m/s^.  I^>max  is 
the  maximum  observed  blade  length. 

Figure  11.  Probability  of  dislodgement,  (P^)  is  plotted  against  water 

acceleration  for  stationary  (left  panels)  and  moveing  (right  panels)  snails 
exposed  to  a  variety  of  water  velocities  (1  to  20  m/s),  (a)  scutulata, 

0.005  m  in  radius;  (b)  T.  emarginata,  0.02  m  in  radius,  (c)  T. 
canaliculata,  0.02  m  in  radius.  Notice  that  no  animal  can  remain  attached 
when  exposed  to  a  wave  of  20  m/s. 

Figure  12.  Predictions  for  the  barnacle,  cariosus,  exposed  to  water 

velocities  of  1  and  20  m/s.  Left  panel:  probability  of  dislodgement  is 
plotted  against  acceleration  for  a  barnacle  0.018  m  in  radius.  Right 
panel:  size-specific  increment  in  risk  is  plotted  against  size  (radius)  for 
a  barnacle  exposed  to  a  water  acceleration  of  500  m/s^.  The  maximum  radius 
observed  is  shown  by  the  vertical  dashed  line. 

Figure  15.  The  Log  of  the  cumulative  probability  of  survival  (Ps>cum)  is  plotted 
against  time  (in  months)  for  a  variety  of  intertidal  organisms.  For  each 
animal,  predictions  are  calculated  using,  the  maximum  observed  size  (L), 
half  that  size  (L/2),  and  twice  that  size  (2L). 

Figure  14.  The  log  of  the  probability  of  survival  is  plotted  against  number  of 
waves  encountered  while  moving  (left  panels)  and  against  time  while 
stationary  (right  panels)  for  three  limpet  species:  (a)  pelta,  (b)  N. 
scutum,  and  (c)  C^_  digitalis. 
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